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[1] DERIVE-Projekte im Unterricht, MU, Jahrgang 41, Heft 4/95 
Friedrich Verlag, 30917 Seelze 

[2] Materialien zum M-Unterricht mit Computer und DERIVE, 

Landesmedienzentrum in Rheinland-Pfalz, 1995, Hofstr. 257, 56077 Koblenz 

[3] Learning Linear Algebra through DERIVE, Brian Denton, 

Prentice Hall, 1995, ISBN 013 122664 9, 351 pages 

[4] Business Calculus today with Spreadsheets and DERIVE, R. L. Richardson, 

Saunders College Publishing, 1996, ISBN 0 03 017554 2, 416 pages 

[5] Resources by Discovery, MAA Notes 27-31, The Mathematical Association of America 

[6] transferts, Cahiers de la Cellule Recherche Innovations Pedagogiques, Numero 7 

Derive: Utilisation d un logiciel de calcul formel 

C.R.I.P. Rectorat de Tacademie de Lille, 20 rue Saint Jacques, 59000 Lille, France 

Dominique Lymer, one of the authors and a DUG member since long sent one copy (170 
pages). It is an interesting collection of work sheets for students with accompanying 
support material for teachers. Merci, Dominique. Dominique has submitted a contribution 
which is one of my favourites for the next issue. 

[7] Matematicas con DERIVE en la economia y la empresa, Alfonso Gonzales, Pareja, 
ra - ma, ISBN 84 7897 201 3 


Exchange for DERIVE Teaching 


Borse fur DERIVE Unterrichts- 

materials in the DNL 


materialien im DNL 

The wheel has not to be invented twice. 


Das Rad muss nicht zweimal erfunden werden. 


I can offer: 

Binomial Theorem, GCD & LCM, System of Coordinates (in English and in Ger- 
man as well), Modelling Word Problems with DERIVE 

DERIVE Days Leeds 

Spring School on Teaching and Learning 
Mathematics with DERIVE and the TI-92 

Trinity and All Saints College, Leeds 
13 - 15 April 1996 


For receiving an announcement please send a self-addressed A4 envelope to: 
John Monaghan, Centre for Studies in Science and Mathematics Education, 
University of Leeds, Leeds LS2 9JT 
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Liebe DUG-Mitglieder, 


Dear DUG Members, 


Gleich zu Beginn mochte ich uns alien zum funf- 
jahrigen Bestehen der DUG gratulieren. Mit dem 
ersten Mitglied aus Namibia , das ich recht herz- 
lich in unserem aller Namen begriifie, sind nun 
alle funf Kontinente in der DUG vertreten. 

Fur das Jahr 1996 steht eine interessante Auswei- 
tung des DNL in Aussicht. Bert Waits von der 
Ohio State University (OSU) und Bernhard Kutz- 
ler von SWHE, zwei TI-92 Spezialisten, haben 
angeregt, eine TI-92 Kolumne in den DNL aufzu- 
nehmen. Sie haben auch versprochen, die ersten 
Beitrage zu liefern und Anfragen zum TI-92 zu 
beantworten. Das „klassische“ DERIVE fur den 
PC wird dabei weiterhin nicht zu kurz kommen. 

Die Beitrage dieser Ausgabe sind schwerpunkt- 
mafiig den Differ entialgleichungen gewidmet. Das 
ist vielleicht fur manche von Ihnen zu einseitige 
Kost. Ich kann Ihnen Abwechslung ankiindigen: 
1996 wird es einen ausgesprochen sportlichen 
DNL geben: eine Projektarbeit ilber das Schi- 
springen, eine Geometrie des Fufiballs, Baseball, 
ein Vergleich von sportlichen Leistungen und 
vielleicht ein Beitrag „ DERIVE und das Tennis- 
netz “ versprechen allerhand. 

Als kleines Weihnachtsgeschenk finden Sie unter 
den Dateien im Unterverzeichnis <STEREO> 
einige Stereogramme, an denen Sie Ihren „magi- 
schen Blick“ testen konnen. Laden Sie die Bilder 
in ein Bildbetrachtungsprogramm (zB 
IRFANVIEW). Sie las sen sich ausdrucken oder am 
Bildschirm betrachten. Ich wiinsche Ihnen viel 
Spafi damit. Die Bilder wurden mit POPOUT- 
LITE erstellt. 

(Geom. Figuren, Geckos + DUG, Geckos, SWHE- 
Logo, DERIVE in Rahmen und Gecko + DERIVE) 

Bitte vergessen Sie nicht, Ihre Mitgliedschaft 
rechtzeitig zu erneuern, wenn Sie der DUG die 
Treue halten wollen. 

Damit bleibt mir nur noch, Ihnen und Ihrer Fami- 
lie in meinem und im Namen meiner Frau alles 
Schone fur das bevorstehende Weihnachtsfest und 
ein gliickliches Neues Jahr 1996 zu wiinschen. 


Herzliche Griifie 



Let me start with a congratulation: It is your 
merit that we can celebrate the 5th year of the 
existence of the DUG. As we have now a 
member from Namibia who joined us two 
months ago we can really say that the DUG is 
represented on all five continents. We will 
express our warm welcome. 

For 1996 I am looking forward to having an 
interesting enrichment in the DNL. Bert Waits 
form the Ohio State University (OSU) and 
Bernhard Kutzler (SWHE) - two TI-92 special- 
ists -suggested including a TI-92 column in 
the DNL. They both have promised to submit 
the first contributions and to answer TI-92 
related requests. 

Some contributions of this issue focus on 
differential equations. Td like to ask our 
friends who are not involved so much in these 
applications for patience. I can promise for 
1996 a really “sportive” DNL: a student’s pro- 
ject on ski jumping, geometry of the soccer 
ball, base ball, a comparison of sportsmen’s 
performances and hopefully an article 
“ DERIVE and the tennis net”. 

Among the files you can find a small Christ- 
mas gift. There are some stereograms in the 
folder <STEREO>. Try your “magic eyes” by 
loading them into any graphics program (e.g. 
IRFANVIEW). You can either print the pic- 
tures or inspect them on the PC-screen. I 
used POPOUT-LITE to create the pictures. 

(Geometric shapes, geckos+DUG, geckos, 
SWHE-Logo, DERIVE in frames and 
gecko+DERIVE). 

Please don’t forget renewing your member- 
ship for 1996. 

My wife and I wish you and your family the 
best for Christmas and a Happy New Year 
1996 

Best regards 



I met Joe Fiedler at a TI-92 workshop in Houston. Unfortunately I didn’t realize that Joe has 
been a DUG member, so I didn’t introduce myself. Joe, you were great. I enjoyed your and 
Wade Ellis' workshop immensely. I had many nice hours at the ICTCM including the DUG- 
meeting. But this workshop was one of the top events for me. Many thanks. 
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The DER1VE-NE WSLE TIER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with a con- 
tents of 30 pages minimum. The goals of 
the D-N-L are to enable the exchange of 
experiences made with DERIVE as well as 
to create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone: 43-(0)2275/8207 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the D-N-L. 
It must be said, though, that non-English 
articles are very welcome nonetheless. 
Your contributions will be edited but not 
assessed. By submitting articles the author 
gives his consent for reprinting it in D-N- 
L. The more contributions you will send, 
the more lively and richer in contents the 
DERIVE Newsletter will be. 


Preview: (Contributions for the next issues): 

Graphic Integration, Probability Theory, Linear Programming, Bohm, A 

LOGO in DERIVE, Lechner, A 

DREIECK.MTH, Wadsack, AUS 

IMP Logo and Misguided Missiles, Sawada, HAWAII 

3D Geometry, Reichel, AUS 

Parallel- and Central Projection, Bohm, AUS 

Vector and Vector Indices Sorting, Biryukov, RUS 

Algebra at A-Level, Goldstein, UK 

Tilgung fremd erregter Schwingungen, Klingen, GER 

Utility for Complex Dynamic Systems, Lechner, A 

Notes on DERIVE 2.6 functions and limits, Speck, NZL 

Ski Jump, a project with students, Scheuermann, GER 

Linear Mappings and Computer Graphics, Kiimmel, GER 

Julia Sets, Kurnmel, GER 

Solving Word Problems with DERIVE, Bohm, AUT 
and 

Setif, FRA; Vermeylen, Belgium; Lymer, FRA; Leinbach, USA, Aue, GER; 
Weth, GER, Wiesenbauer, AUT; Keunecke, GER, Weller, GER, ... 

and messages from the derive-news@mailbase.ac.uk 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag.Josef Bohm 
Herstellung: Selbstverlag 
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across a concept (or topic) which can be investigated from an educational point of view, and 
which has, in their opinion the potential to contribute to our understanding of how students 
learn mathematics via a computer. 

Regards. Martin Lindsay, e-mail: Martin Lindsay@vut.edu.au; FAX: ++61 36884050 

DNL: I sent some own materials to Martin and I hope that there are some other DUG members who 
will support his work. Are there any suggestions for projects? It could be interesting to compare pro- 
jects done in different countries und under various circumstances and conditions. Hopely we will hear 
about results. Additionally I can recommend the Resources for Calculus Collection from the MAA, 
which I bought in Houston last week and transferts (in French). You will find more information in the 
Book Shelf. 

In the last DNL you could find an interesting request about a square root simplification. I received two 
respective answers, one from the Fachhochschule Osnabriick and another one from Albert Rich, 
SWHH, which explains the issue. 

FHS Osnabriick 

The fact described by Mr. Propper in DNL#19 has more odd aspects. If one having defined 
the interval [-1,1) for x, tries to simplify expression #1 he will obtain the expected result (#5). 
Only expression #2 (derived by DERIVE) will not simplify. Although DERIVE recognizes #1 
and #3 as equal, they are treated differently when they are simplified. 

The next example shows the same: the declaration of the interval -1 < x 1 is considered in 
simplifying expressions #1 , #9 and #1 1 . Simplifying #1 3 (included #2) fails again. 


#1: 

(l + x) 

7((1 + x).(l - x)W 

l 1 - X J 

User 

#2: 

x + 1 'I 

7 -7((x + D-d - x)) 

U - X J 

Si mp(#l) 

#3: 

2 (l + x) 

7(i - x ) -7 

U - x J 

Expd(#l') 

#4: 

x :e Real [-1, 1) 

User 

#5: 

x + 1 

Simp(#l) 

# 6 : 

2 

7(1 - x )-7(x + 1) 

Simp(#3) 


7(1 - x) 


F(x) := 

If x > 2 

#7: "false" User 

"true" 

"unknown" 


#8: F(x) := true 


Simp(#7) 
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PART I T I ON S_AUX ( n , m , mn ) := 

IF (m=l , 1, SUM ( PARTI TIONS_AUX (n-k_,m-l, k_) , k_,mn, FLOOR (n,m) ) ) 
PARTITIONS (n) : =SUM (PARTITIONS_AUX (n, m, 1) ,m, l,n) 

PARTITIONS (4) =4 

Unfortunately, for n greater than about 50 the above functions take too long because 
of explosive fan-out on the recursive calls to PARTITIONS_AUX: Are there any DUG mem- 
bers who can come up with an efficient definition for PARTITIONS that Soft Warehouse can 
include in NUMBER. MTH? Aloha. 

(Are there? Hannes and others, this is a real challenge! I have put these functions and the next ones 
together in one file NUMBEXT.MTH. You can merge this file to NUMBER.MTH and so improve 
your NUMBER.MTH. Josef) 

Albert Rich #2 

Dear Josef, 

Enclosed is an article that appeared in a recent Science News about consecutive 
prime numbers that you might find interesting. I thought that these seven 97-digit primes 
would make a good test of DERIVE's NEXT_PRIME function. The enclosed printout shows 
that DERIVE does correctly recognize primes of this size. 

Since my October 16, 1995 letter we have simplified and renamed the function for 
computing the partitions of a number. This is how my comment in this letter should read: 

PARTS(n) simplifies to the number of decompositions of n into integer summands 
without regard to order. For example, 4 = 1+ 3 = 2 + 2 = 1 + 1+2 = 1 + 1 + 1+1so 
PARTS (4) = 5. 


PARTS_AUX (n, m) -IF (n<2m, 1, 1+2 (PARTS_AUX (n-k_, k_) , k_,m, FLOOR (n, 2) ) ) 
PARTS (n) —IF (n<l , 0 , PARTS_AUX (n, 1) ) 

PARTS (4) =5 
PARTS (50) =204226 

(needs 8.05 sec) 


Recent Derive versions have another PARTS-function implemented. It looks very bulky but it 

is much more efficient than the earlier one - and it works without the auxiliary function: 

PARTS (n) :=IF (n<2, 1, FLOOR (APPROX (2 (-1 /tW (k_/2) 2 (IF (GCD (h_, k_) =1, 

COS (tt( 2 ( (i_/k_-l/2 ) (MOD (i_h_/k_) -1/2) , i_, 1, k_-l) -2nh_/k_) ) , 0) , 
h_, 1, k_) 2V6e A (- tW (24n-l) / (6k_) ) (e A (tW (24n-l) / (3k_) ) 

(6k_-7iV (24n-l) ) -6k_-/W (24n-l) ) / (k_(24n-l) A (3/2) ) , 

k_, 1, Vn/LOG(n, 11) ) ,LOG(l/ (4nV3)e A (tW (2n/3) ) , 10) +5) +1/2) ) 

PARTS (50) =204226 

(needs 0.000 sec!! Josef) 
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DISTINCT_PARTS(n) simplifies to the number of decompositions of n into distinct in- 
teger summands without regard to order. 

For example, 4 = 1 + 3 so DISTINCT_PARTS(4) = 2. 

DISTINCT_PARTS_AUX (n,m) -IF (n<2 -m, 1, l+£ (DISTINCT_PARTS_AUX (n- 
k_, k_+l ) , k_, m, FLOOR (n, 2) ) ) 

DISTINCT_PARTS (n) —IF (n<l , 0 , DISTINCT_PARTS_AUX (n-1, 1) ) 

DISTINCT_PARTS (100) = 444793 

(needs 16.5 sec) 

The function which is implemented now works also without any auxiliary function. Compare 
the calculation time. 

DISTINCT_PARTS (n) -£ (IF (n-1 6i_+l=FLOOR (V (8n-16i_+l) ) / '2, PARTS (i_) , 0) , 

i_, 0, n/2 ) 

DISTINCT_PARTS (100) =444793 

(needs 0.047 sec) 

All for now, 

Aloha, Albert D. Rich, Applied Logician 

(Albert's letter closes again with his challenge to improve these functions. Find now the included file 
for Harvey Dubner's first of seven consecutive primes separated by 210. That is a new record, the 
previous had been numerous examples involving six consecutive primes in arithmetic progression. 
Albert refers to an article from SCIENCE NEWS, vol. 148, September 1995. The article starts 
"Searches for patterns among prime numbers have long served as stiff tests of the ingenuity and per- 
severance of mathematicians." You can obviously see that not only our column writer Johannes 
Wiesenbauer is dealing with prime numbers. H. Dubner and H. Nelson ended up using seven com- 
puters, running continuously for about 2 weeks, to find the sequence. Now they are thinking about 
going to eight consecutive primes. They estimate that it would take 20 times longer - at least 2.5 com- 
puter years - to accomplish this search on their souped-up personal computers. Josef) 

"Harvey Dubner's first of 7 consecutive primes separated by 210:" 

p: =108 95334 3124 7 05 931 08757 8037 8 922 957 732 90 803 64 92 9 931 381 953852 
13105561742150447308967213141717486151 

"The following verifies the difference of the 7 consecutive primes" 

"and that DERIVE ' s primality is working good:" 

v : =ITERATES (NEXT_PRIME (n) ,n,p, 8) 

"You might simplify #5 to see the 8 huge prime numbers" 

VECTOR (v SUB (n +l)-v SUB n , n, DIMENSION (v) -1 ) 


[210,210,210,210,210,210,120,52] 


p 
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Playing-Cards Shuffling with DERIVE 

Benno Grabinger, Neustadt, Germany 


In DNL#12 Mr Chuan from Taiwan posed the following question: 

Is it a surprise that after perfectly shuffling only 8 times a deck of 52 cards, the origi- 
nal position of the cards is restored? 


With the help of DERIVE it is no problem to understand what is going on. The perfect shuffling is 
described by the permutation 

_fl 2 3 4 ... 49 50 51 52^ 

P ~{\ 27 2 28 ... 25 51 26 52/ 

It is well known that each permutation is a product of cycles (having no common elements). The func- 
tion CYCLE (p) creates a product of cycles which represents the permutation p\ 


Applying this function to the given permutation p leads to: 


[[1], [27, 14, 33, 17, 9, 5, 3, 2], [28, 40, 46, 49, 25, 13, 7, 4], [29, 15, 8, 30, 41, 21, 11, 6], [31, 16, 34, 
43, 22, 37, 19, 10], [32, 42, 47, 24, 38, 45, 23, 12], [35, 18], [36, 44, 48, 50, 51, 26, 39, 20], [52]] 


It can easily be seen that this product consists of 6 cycles of length 8 and one cycle of length 2. There- 
fore, if p has been applied 8 times the original distribution will be restored. 


#1: 

■CYCLE. MTH by 

Benno Grabinger, 

199 r >H 



#2: 

n := 52 





#3: 

u := VECTOR (k, 

k, 1, n) 




#4: 

V := VECTOR ^IF 

|mod(i, 2) = 0, 

i n 

i + 1 > 

. i, 1, nj 

— + — p 

2 2 

2 > 


#5: p := [u, v] 


1 

2 

3 

4 5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

1 

27 

2 28 3 

29 

4 

30 

5 

31 

6 

32 


7 

33 

8 

34 

9 

35 

10 

36 

11 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


33 

34 

35 

36 

37 

38 

39 

40 

41 

37 

12 

38 

13 

39 

14 

40 

15 

41 

16 

42 


17 

43 

18 

44 

19 

45 

20 

46 

21 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 











47 

22 

48 

23 

49 

24 

50 

25 

51 

26 

52 












#7: VALUE(i, p) := ELEMENT(ELEMENT(p , 2), i) 

#8: Z(S, p) :i DELETE_ELEMENT(ITERATES(VALUE(], p), ], s), 1) 

EQUAL(V, w) := 

If DIMENSION(w) = 0 
0 

#9: If ELEMENT(v, 1) = ELEMENT(w, 1) 

1 

EQUAL(V, DELETE_ELEMENT(w, 1)) 

IN(v, z) := 

If DIMENSION(v) = 0 
0 

#10: If EQUAL( ELEMENT (v , 1), z) = 1 

1 

IN(DELETE_ELEMENT(v, 1), z) 
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CREATE(p , i, v) := 

If i r II + 1 

v 

#11: If IN(v , Z(i , p)) = 1 

CREATE(p, i + 1, v) 

CREATE(p, i + 1, APPEND(v , [Z(i , p)])) 

#12: CYCLE(p) :: CREATE(p , 1, []) 

#13: CYCLE(p) 

#14: [[1], [27, 14, 33, 17, 9, 5, 3, 2], [28, 40, 46, 49, 25, 13, 7, 4], [29, 15, 8, 

30, 41, 21, 11, 6], [31, 16, 34, 43, 22, 37, 19, 10], [32, 42, 47, 24, 38, 45, 
23, 12], [35, 18], [36, 44, 48, 50, 51, 26, 39, 20], [52]] 


Two nice graphics from Belgium 



G P Speck, the author of the contribution “Mueller's Method” lives in Wanganui, New Zealand. 
The Wanganui River is the main river on the Northern Island of New Zealand. One of the famous 
places in the Wanganui region is the “Bridge to Nowhere”. Regards to wonderful New Zealand 
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SOME IMPROVEMENTS ON THE RESOLUTION OF 
ORDINARY DIFFERENTIAL EQUATIONS 

Jose-Luis Llorens Fuster 
Universidad Politecnica de Valencia 
Departamento de Matematica Aplicada 
46071 -VALENCIA, SPAIN 


First order equations 

The general solution of the first order differential equation 

p(x, y) dx + q(x, y) dy = 0 [a] 

can be obtained using function DSOLVEl_GEN(p,q,x,y,c), which is incorporated in the 
ODE1.MTH file (see [1], chap. 5, p.333; [2] chap. 9.6, p.252), always that this equation be: 

• Separable 

• Linear 

• Homogeneous 

• Exact 

• Equation having an integrating factor which depends only on x or only on y. 

When the differential equation is not one of the previous type, if we simplify the men- 
tioned function we obtain “inapplicable” . If we want to obtain the particular solution of this 
differential equation satisfying the initial condition (xo,yo)> we can use the 
DSOLVEl(p,q,x,y,xo,yo) function. 

Thus, the general solution of the equation (taking from [1], p. 353, ex. 3) 


Ex. 1 : y' sin 2 x = cos 2 y 

can be obtained simplifying the expression DSOLVEl_GEN(cos 2 y, -sin 2 x). The general solu- 
tion will be obtained depending on the constant ceR. 

In the ODE1.MTH file (which is now FirstOrderODEs.mth, Josef) there are some in- 
dependent functions for these types of differential equations: 

SEP ARABLE_GEN (p, q, x, y, c) for the separable differential equation having the 
form y' = p(x) q(y). 

LINEARl_GEN(p, q, x, y, c) for the linear DE having the form y' + p(x)y = q(x). 
HOMOGENEOUS_GEN(r, x, y, c) for the homogeneous, 

EXACT_GEN(p, q, x, y, c) for the exact, 
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INTEGRATING_FACTOR_GEN(p, q, x, y, c) for these equations having an inte- 
grating factor which depends only on x or only on y. In the last three functions we suppose 
that the differential equation is written in the implicit form [a]. The DSOLVEl_GEN func- 
tion chains the previous functions using the corresponding tests to identify each type of the 
equation. For example, if the differential equation satisfies: 

dp _ dq 
dy dx 

then it is an exact one. Thus the definition of the function contains a sentence of the type 
IF(DIF(p,y) = DIF(q,x), EXACT_GEN(p,q,x,y,c),...). From this we can conclude that the 
efficiency of this function depends not only on the previous five definitions but also on the 
corresponding test. 

On the other hand, it is possible to modify the way how to present the differential 
equation [a]. For example, we can write it as: 

dy _ 

dx q(x,y ) 

So, in order to solve it we simplify the expression 

DSOLVE l_GEN(-p(x,y)/q(x,y),-l ) 

and obviously we expect to obtain the same solution (or an equivalent one). However, the 
following illustration shows the behaviour of DERIVE (v. 3.01) with respect to the separable 
differential equation of example 1 : 



Obviously, the bug is not in the resolution of the equation (because, as you can see, 
the SEP ARAB LEGEN function works correctly) but in the test. But this is not an excep- 
tional example: the following differential equations (taking from [1], p. 353, ex. 15, 18, 43, 45 
and 46) are also of separable form: 

Ex. 2: x 3 y' = y 2 (x - 4) 

Ex. 3: xy'(2y-l) = y(l-x) 

Ex. 4: sec 2 xtan y dy + sec 2 y tanx dx = 0 

Ex. 5: x 2 yy r = (l + x)cscy 

Ex. 6: sinj;cos 2 j;dx + cos 2 xdy = 0 
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Z 3 

DSQLUE 1 JBEN (y ■ Cx - 4), - x > = "Inapplicable" 

DS0LUEl_6ENty <1 - xJ, - tZy - 1> x> = "inapplicable" 

Z Z 

DSOLU E1_GEN ( SIC C y ) -TrtMCxl, SECtxl -TOHtyJJ = “inapplicable" 


BSDLUEl_GEM C ( 1 + x? CSCtij), x yj = "inapplicable’' 

g E 

DSQLUEl_6EM(SIH(y) COS£y> , COS(x> 1 - "inapplicable' 


Before following Llorens Fuster's very valuable suggestions from 1995 let's look at 
DERIVE 's solutions from today, Josef 


The way how DERIVE 6 solves these differential equations: 


2 3 

#1: DSOLVEl_GEN(y • (x - 4) , - x ) 


1 

y 


2 ^ 
c-x - x + 2 


#2: DSOLVEl_GEN(y ■ (1 - x), - (2-y - l)-x) = (LN(y) - 2-y = - LN(x) + x - c) 

2 2 2 2 
#3: DSO LVE1_G E N ( S EC (y ) -TAN(x), SEC(x) -TAN(y)) = (SIN(y) = 2-c - SINGO ) 


#4: DSO LVE1_G E N ( (1 + x).CSC(y), x -y) = 


ox + 1 


yCOS(y) - SIN(y) = LN(x) 


2 2 

#5: DSOLVEl_GEN(SIN(y).COS(y) , COS(x) ) = 


r 

r 

f y T 

\ i > 

- LN 

TAN 


= TANOO c 



^ 2 

J COS(y) J 


The TI-92, Voyage 200 and TI-NSpireCAS are performing pretty the same: 



We continue with Llorens Fuster from 1995: 
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The way to solve this problem is suggested from the first example. In all these cases, 
DERIVE seems to recognize more easily that the differential equation is separable if it is writ- 
ten in the explicit form [b] than in the implicit one [a]. So, we propose the following modifi- 
cation in the ODE1.MTH file: 


a) We change the name of the functions DSOLVE1GEN and DSOLVE1 respectively by 
SOLVE 1_GEN and SOLVE 1. 

b) After that we add the following functions: 

DSOLVEl(p, q, x, y, xO, yO, a_) : = IF("inapplicable" = a_ : = SOLVEl(p, q, x, y, xO, yO), 
SOLVE 1(- p/q, -1, x, y, xO, yO), a_, a_) 

DSOLVEl_GEN(p, q, x, y, c, a_) : = IF("inapplicable" = a_ : = SOLVE l_GEN(p, q, x, y, c), 
SOLVE 1_GEN(- p/q, -1, x, y, c), a_, a_) 


The meaning of these new functions is obvious: If the original function DSOLVE1 
(now named SOLVE 1) gives “ inapplicable ” then we try to apply it again writing the differen- 
tial equation in the form [b]. Thus we solve all the examples from above: 


2 Z 

JSOLUEl_«EM<CQS(y) , - SIHtx) ) 

- COKx) - TfiN(y) = c 

2 3 

DSOLUEljGENCy -(* - 4), - x > 

12 1 

x 2 y 

X 

DSOLUEl_GEHty Cl - x>, - tZy - !)-■> 
LHfx-y) - x * 2-y “ c 


2 Z 

®SOLUEl_GEMCSECCy > TfiNf*}, SECtxl TAMlyH 

z z 

SIMfx) SIMfyi 


z z 

z 

DSOLUEi_GEMai + x) CSCCyJ, x 

1 

- LMCx) * y COS(y) * SIMy) * — - c 

x 

Z Z 

IVSOLMEl.GEmSIHCy) COSCy) , COStxJ ) 

- yiUflll - 1WCx> - 1 — = c 

1 1 Z JJ COSCy ) 


We have suggested to Soft Warehouse to include “officially” this modification in the 
ODE1.MTH file, although some strange examples are still remaining: 

Ex. 7: (xy 1 - y 2 + x - l)dx + (x 2 y - 2x y + x 2 + 2y - 2x + 2 ) dy = 0 

(taking from [1], p. 354, ex. 31) which admits the integrating factor //( v) = e ~ 2arctanv . How- 
ever, in this way the general solution is not obtained. The new definition of the DSOLVE1 
function does not act because the given result is not “inapplicable” . DERIVE identifies cor- 
rectly the type of differential equation but is not able to obtain the solution because it cannot 
evaluate the integrals involved: 
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Z 2 Z 2 

BSni.UE1_GEN( it -y - y * * - 1, * y - 2 x y * * * Z y - Z x + 2) 


Z Z 1 

£ - ATAHCy > x ■ <y + 11 Z | f ZATMMy) f 2- 

* - 1 — - — - — — - x + 1)1 * 2 J y dy + Z-J £ 

L Z * 


ATAfl 


D50lUEI_6En 


2 2 

x-y - y ■+ x - 1 


Z Z 

* y-Z x y + Z-y - Z-x + z 


-* -1 


z z 

LNCIm - Z x * Z) Cy * m 
- ATftH (y ) - ^ = c 


As it is shown in the third expression in the previous illustration, it is solved if we 
write the equation in form [b]: This is due to the fact that now the integrating factor is 

/+1 


M = ' 


T + l 


, which is leading to an exact differential without appearing the mentioned inte- 


grals. Obviously this example is very rare. To solve this problem we need to modify more 
deeply the ODE1.MTH file. 


As you can see even DERIVE 6 is not perfect in solving ODEs. It needs the same rewriting 
as in 1995. 


2 2 2 2 

#1: DSOLVELGEN(x -y - y +x-l, x-y - 2-x-y + x + 2-y - 2-x + 2) 


#2: c 


2 -ATAN(y) 


2 2 

x -(y + 1) 2 

x-(y + 1) 


rr 2-y 
e 


+ 2-SUBST 


COS(y) 


dy, y, ATAN(y) 


2-SUBST 


2-y 

G -SIN(y) 


dy, y, ATAN(y) 


COS(y) 


= c 


#3: DSOLVELGEN 


2 2 

x-y - y + x - 1 


2 2 

x -y - 2-x-y + x + 2-y - 2-x + 2 


-1 


2 2 
LN((x - 2-x + 2) • (y + 1)) 

#4: ATAN(y) + = -c 

2 


It is interesting that the Tl-family shows no 
problem to solve this “rare” kind of differen- 
tial equation without hesitating. 


ifl^T FF? Y F]:t T FHt y ft Y Fe^ Y '| 

I t i-— |fl 1 gebra |Ca 1 c lother |Pr gm 1 0 |C 1 ean Up | | 

1 1 1 1 1 — -i 1 1 — _ -i- _ i i 1 — 1 


lM-y" - y*- + x - ij 
2-x-y + x 2 + 2- y- 2V 


■ deSoloe y 1 = 


ln(y 2 + l) 


2 

x ■ y ■ 


+ tan-^y) = £6 - 


ln(x 2 - 2-x + 2) 


deSoluefy 1 (x 2 y-2xy + x 2 + 2y-2x> 
ln(y 2 


^ + tan-i(y) = S3 - ln[ * 2 ^ * + 2) 


... -2x+2 >= ~ <x*v A 2-<j A 2+x-l > ■ x ■ v>l 

MAIN RAD AUTO FI1NC ZW 
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Finally we show an example that is neither a limitation of this function nor in general 
of the program. The linear differential equation 

Ex. 8: y' + ycosx = sin2x 

([1], p. 366, ex. 56) is not correctly solved if we don't select Manage-Trigonometry- 
Expand: 

: DSOLIIE 1_GEN (S IN (2 - x ) - y-COS(x), -1} 

SlMtxJ SIM(x) 

t SlfltZ x) dx - y 6 = c 

Tr IgnnnMfttry := Expand 

31111x1 

CZ SIWCx] - y - ZJ = c 

This is no problem for DERIVE 6. 

But it is a problem for the Voyage 200 - and 
TI-NspireCAS as well - which can be re- 
solved expanding the trig expression before 
integrating the differential equation. 


S3 


Y F]:t T FHt Y ft Y FT^ Y i 
fllgebra|Calc|Other PrgniIO|Clean Up| || 
fla^y^fe--'^ 2 " " — 

1 deSolue(y 1 + y ■ cos(x) = sin( 2 ■ x) , x , y) 
y = *■ ^ + (29-«> 

'tExpancKy' + y ■ cos(x) = sin(2 ■ x)) 

cos( x) ■ y + y 1 = 2 ■ s i n( x) ■ cost. x) 
1 deSolye(cosCx) ■ y + y 1 =2 ■ sin(x) ■ cos(x) ? >► 
y = eiQ-f + 2(sin(x) - 1) 


desolue<ans<l>,x,y> 


#5: 


# 6 : 


DS0LVE1_GEN(SIN(2 • x) - y.COS(x), -1) 


SIN(x) 

e • (2 * SIN(x) - y - 2) 



This is the recent form of DSOLVE1_GEN how it is implemented in DERIVE 6. Jose-Luis 
suggestion for improving the version from DERIVE 3 was obviously accepted by Soft Ware- 
house and extended to include other cases, too. 


DSOLVELGEN(p, q p x, y p c p a_) := 

If "inapplicable" = a_ :z INTEGRATING_FACTOR_GEN(p p 
If "inapplicable" = a_ :z HOHOGEN E 0 US_GEN ( - p/q p 
If "inapplicable" = a_ := SEP_GEN(p, q p x p y p 
GEN_HOH_GEN(- p/q, x p y p c) 
a_ 
a_ 
a_ 
a_ 


q, x p 
x, y, 

c) 


y. c) 
O 


a_ 

a_ 
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Second order equations 


When a second order differential equation y" = u(x,y,y') is incomplete, i.e., when 
some of its terms in x, y, or y ’ do not appear, sometimes it is possible to solve it directly or 
reduce it to two first order equations. So: 


a) 


y and x do not appear: the equations are of the form y" = / (y'). The parametric equa- 
tions of the general solution as a function of the parameter v and the constants c and k , 
are: 


x = 


dv 

W) 


+ c, y = 


vdv 

W) 


+ k. 


b) y’ and x do not appear: the equations are of the form y" = f (y). The general solution 

can be written in one of the following expressions as a function of the constants c and k : 
r 

dv 


x = 


J2J/00# 


+ k 


r + c, y 


= J ^2j/(y)dy+I dx + c. 


c) v does not appear: using the change v = y' we transform it to a first order differential 
equation (in v\ v, x). The general solution of this equation is, obviously, another first 
order equation. 

d) x does not appear: the changes y' = v, y” = v' ■ v, transform it to a first order differential 
equation with the variables v (in place ofy) andy (in place of x). The general solution of 
this equation is, again, another first order equation. 

It is well known, that in the ODE2.MTH file are only two functions to solve the dif- 
ferential equations presented in a) and d). In particular: 


AUTONOMOUS(r,v). We apply it to solve equation d), written in the form 
y" = Hy, v) with y ' = v. DERIVE simplifies it solving for v’ after applying the previous 
change. 


AUTONOMOUS_CONSERVATIVE(r,x,y,xo,yo,vo). We apply it to solve equation 
d), written in the form y" = r(y), where the initial conditions y(x 0 ) = y 0 , y'(x 0 ) = v 0 are sup- 
posed. So, it cannot be used to obtain the general solution of the equation. 


Our objective is to define auxiliary functions to obtain the general solution of the dif- 
ferential equations of kinds a), b), c) and d). 


The following shows what Jose-Luis proposed in 1995. Then we will have a look how the 
utility file SecondOrderODES.mth is supporting solving second order DEs in our times. 
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3 For case a), i.e. y" = u(y ' ) = u(v), it is very easy to define the function: 

FALTAYX (u,v,x,y,c,k):= [x=INT (1/u, v) +c, y=INT (v/u, v) +k] 


3 Similarly, for the incomplete equations of kind b), i.e. y” = u (y ) : 

FALTAXV (u , v) == ["Selec. : ",x=INT ( (2f (u, y) +k) A (-1/2 ) , y) +c, 
y=INT (2f (u, y) +k A (1/2) ,x)+c] 


In this way we obtain two equivalent expressions of the general solution. It is possible 
that one of the two solutions leads to an unsimplificable integral, depending on the form of 
the equation; for these reasons we present the two expressions. 

3 Finally, for the equations of type c) and d), i.e. y" = u(x,y) or y" = u(y,v) respec- 
tively, we have to define in the first step two auxiliary functions with the following 
purpose: 

to solve the first order differential equation obtained after y’ = v as we stated 
previously, 

to solve for v in the general solution of this differential equation. 


AU1 (u,v,x,y,k) - SOLUTIONS (DS0LVE1_GEN (u,-l,x,v,k),v) 


AU2 (u,v,x,y,k) -SOLUTIONS (DS0LVE1_GEN (u/v, -1, y, v, k) , v) 

(I had to adapt the version from DERIVE 3 to DERIVE 6 because of another form of output 
for the SOLVE command.) 

Now we define the corresponding main functions in which we apply the 
DSOLVE1GEN function in the obtained solution with each auxiliary function. As it is pos- 
sible to obtain two or more different expressions when we solve for v, we use the VECTOR 
command: 


FALTAY (u,v,x,y,k,c) : = 

VECTOR (DSOLVEl_GEN( (AU1 (u,v,x,y,k) ) |i,-l,x,y,c) ,i, 
DIM (AU1 (u, v, x, y, k) ) ) 


FALTAX (u , v,x,y , k , c) : = 

VECTOR (DS0LVE1_GEN( (AU2 (u,v,x,y,k) ) |i,-l,x,y,c) ,i, 
DIM(AU2(u, v, x, y, k) ) ) 

The efficiency of these functions depends on the following three factors: 

- The fact that the first order differential equation obtained after doing the substitution 
of the variables belongs to these kind of equations which are recognized by DSOLVE1GEN; 

- The possibility of solving for v in the solution of this equation; 

- The fact that the new differential equation obtained in this process is again of a type 
of equation which is recognized by DSOLVE1GEN. 
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In the following functions we study the way how DERIVE identifies the kind of in- 
complete equation and, as a consequence, applies the right function to solve it. 

TESTX (u , v) :=IF ( 3 (u, x) =0 , FALTAX (u, v) , "inapplicable" , 
"inapplicable" ) 

TESTY (u,v) -IF (3 (u, y) =0, FALTAY (u, v) , TESTX (u, v) , TESTX (u, v) ) 


TESTV (u,v) —IF (3 (u, v) =0 AND 3 (u, x) =0 , FALTAXV (u, v) , 
TESTY (u,v) , TESTY (u, v) ) 

ODE2I (u,v) —IF (3 (u, x) =0 AND 3 (u, y) =0 , FALTAYX (u) , 
TESTV (u, v) , TESTV (u, v) ) 


DSOLVE2I (u,v,w) ==ODE2l ( (SOLUTIONS (u,w) ) |1, v) 


We apply the ODE2I function for the second order equation y" = u(x,y,v). The result 
is “inapplicable” if the equation does not belong to any of the types a) , b) c) or d). In other 
case we obtain the solution with the restrictions noted previously for the kinds c) and d). 


We apply the DSOLVE2I function for the second order equation u(w,v,x,y) = 0, 
where w = y”. Thus, is the main function, where it is not necessary solving w = y” before. 
Given an equation, it is sufficient to substitute w = y", v = y , as one can see in the following 
examples (taking from [1], p. 357): 


Ex. 9 (type c): 


y 

y = xsmx 
x 





2 


" 1 


k-x 

DSOLVE2I 

w - -v = x-SIN(x) 

= 

COS(x) + x-SIN(x) + + y = -c 




2 


Ex. 10 (type c): y ' = y"x + y" 2 


2 

DSOLVE2I(v = wx + w ) 





- 2-k -k 


y = 

IF 

x > 0, IF 

k k e e |x| 

(x-e < -1 a x < 0) v (x-e < 1 a x > 0) , 

L 4 2 1 

dx + c 

f 


' ( -2-k -k 
e e -x 


L 4 2 


2 -k 

X ■£ 

dx + c = - 

4 


- 2-k 

x-e 

+ + c 

4 


Jose-Luis Fuster presents another result. So 
I'd like to double check the solution. I define 
a function y_(x) and substitute in example 


10. 


- 2-k 2-k 

x-e x -e 

y_(x) := + c 

4 4 


2 

y _ ' (x) - (y _ ’ r (x) ■ x + y_' r (x) ) = 0 


This solution seems to be correct. 


Ex. 11 (type c): 


//=(i +/ 2 r 


2 1/2 

DS0LVE2I(v ■ w = (1 + v ) ) 


2 2 

2 LN(^(v + 1) + v) v-^(v + 1) 

x = c + ^(v + 1) , y = - + k + 


Ex. 12 (type a): 


y'-y' 2 = 1 


DS0LVE2I(w - v = 1) = 


LN(v + 1) 

x = ATAN(v) + c, y = + k 


Ex. 13 (type c): 


(\-x 2 )y"-xy'-2 = 0 


DS0LVE2I((1 - x )-w - x-v - 2 = 0) 


f 2 2 

U-NC/(x - 1) + x) + 


k.LNGftx - 1) 


+ x) + y = -c] 


Ex. 14 (type a): 


y 




DSOLVE2I(w + ^(1 - v )) 


2 T 2 

= Lx = c 


ASIN(v) , y = 4(1 - v ) + k. 


Ex. 15 (type d): 


yy" = y ' 3 


DSOLVE2I(y-iif = v ) = [y-LN(y) + x - y- (k + 1) = -c] 


Ex. 16 (type d): 


y' 2 =yy' 


2 f 

-k 1 

DSOLVE2I(v = y-w) = L 

x-e - LN(y) = cj 

Ex. 17 (type d): 

y'-y ,3 +yy» = 



3 


3 3 

y 2 


DSOLVE2I(w - v + y-v ) = 

2-x - + y - 2-k-y = c 

3 



Ex. 18 (type d): 


y' = y“lny 


DSOLVE2I(w = v ■ LN (y) ) = 


2 2 
y ■ LN (y) 3-y 

+ x - k-y 
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Ex. 19 (type a): 2 y'y" = 1 + y' 2 



Ex. 21 (type b): y" = 9 y 



As the illustration shows, the differential equation in the next example can also be 
solved taking into account that it is a linear equation. The next example corresponds to a not- 
incomplete linear equation (although y ’ does not appear): 

Ex. 22: y”-y = cos2x-2sin2x 

DS0LVE2I(w - y = C0S(2-x) - 2-SIN(2-x)) = inapplicable 

The two final examples show some limitations. In the first, the expression of the ob- 
tained solution seems not to be very satisfying. In the second (taking from [2], p. 261) we do 
not obtain a solution. In both cases the problem is caused by the resolution of the first order 
equations, i.e. from DSOLVE1GEN. 


Ex. 23 (type d) : 1 + y ,2 =yy" 











D-N-L#20 


Comments on the Differential Equations 


p21 


Ex. 24 (type d): y y" + 1 + y' 2 = 0 
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Some comments for the revised version 

I'd like to show how DERIVE 6 and the TIs are treating the ODEs of 2 nd order. First of all I'll 
solve one or the other of the presented differential equations in the traditional way. This 
might support the students' understanding of the implemented CAS procedures. 

Let's start with Ex. 9 (type c): 

(1) y"-—y' = x sinx 

x 

We substitute y’ = v and y” = V in order to reduce the given DE to a linear DE: 

(2) v'- — v = xsinx. 

x 

The standard technique for solving linear DEs of the form / + p(x)y = q{x) is to find an inte- 
grating factor p(x) and then multiply both sides of the equation by this factor: 


(3) 


X 

— • v' — y • v = sin x 
x x 


J--U 

p(x) = e 1 x =e 



= sinx 


| Integrate wrt x 


— -v = -cosx + A: 
x 

v = -xcos x + kx — > y' 


, . kx 2 

-xcosx + kx — > y = -xsinx-cosxH he 

2 


The last step is easy work: We can integrate directly applying integration by parts. 
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I proceed with Ex. 12 (type a): 


y' = y' I 2 + 1 


found in my (old) textbooks that substituting \/for y’ might be successful. So let me try: 


, d y dy 

dx v 

, „ dv v-dv 

then y = — = 

dx d y 


I perform the substitution in (1) and integrate the 1 st order DE by separation of variables: 


Integrate 


vdv . vdv 

— = v +1 -> dy = — - 

dy v +1 

ln(v 2 + 1) . 

y = — - + k 

2 


dv , dv 

— = 1 + v 2 — » dx = t- | Integrate 


x = tan v + c 

This is the parameter representation which is given as result of example 9. Finally we can try 
to have an explicit form of the solution: 

From (4): v = tan(x-c) and substituting in (3): 

1 ')+*=— ln|cos(jc—c)| + * . 

2 2 l^cos (x-c) J 

But I can also try solving the DE directly by integrating twice: 

dv o dv 

— = 1 + v 2 — >(bc = ^ | Integrate 

dx 1 + v 

x + c = tan -1 (v) v = y f = tan(x + c ) 

— = tan(x + c) | Integrate 
dx 

y = -ln|cos(x + c)| + ^ 


Let's have a type b) example: Ex. 21: (I don't recognize the DE as a differential equation 
with constant coefficients and apply again an appropriate substitution for reducing the order 
of the equation. 


, i „ dv v-dv 

v = v and v = — = 

dx dy 

9ydy = vdv \ Integrate 


<lv py-t, 

lnl\/9v +k. + 3v) 


= dx | Integrate 


9 y 2 +k l =v 2 v = J9y 2 + k 2 


I am trying to bring the solution in its explicit 

form. 
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I would leave it to the students to compare this solution with the solution given on page 19 as 
the result of Llorens Fuster's function. You can also double check the solution by substituting 
into the given DE. 

I must say that I enjoyed comparing the traditional way resolving the equation(s) and their 
outcomes with the CAS results. 


How does DERIVE 6 perform? 

DSOLVE2 (p, q, r , x, cl , c2 ) simplifies to an explicit gen- 
eral solution of the linear second order ordinary dif- 
ferential equation 

y" + p (x) -y 1 + q(x) *y = r (x) 

in terms of arbitrary constants cl and c2 . Note that 
the last two arguments can be omitted if they are vari- 
ables and you are satisfied with the names cl and c2 . 

If no method applies or the equation cannot be con- 
verted to an equivalent one having a p and a q that are 
independent of x, DS0LVE2 returns the word "inapplica- 
ble”. (Online Help) 
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Example 9: 


#1 : DSOLVE2 


SOLVE2^ , 0, x-SIN(x)j = inapplicable 

SOLVEl_GEH^ v - x-SIN(x), 1, x, v j = 

OLVE^COSCx) + = c, vj = (v = ox - x-COf 


#2: DS0LVE1_GEN 


= COS(x) 


#3: SOLVE COSCx) + = c, v = (v = c-x - x-COS(x)) 


#4: DSOLVEl_GEM(c ■ x - x COS(x), -1, x, y, k) = |CQS(x) + x-SIN(x) + y = -k 

As you can see, DSOLVE2 does not apply, one has to perform the reduction of the order. 

AUTONOMOUS (r, v) simplifies an expression for dv/dy, 
given an autonomous equation y" = r (y, v) with v repre- 
senting y ' . 


Example 12: 

2 

2 v + 1 

#5: AUTOHOHOUS(v + 1) = 


#G: DS0LVE1_GEH , -1, v, y = 

2 

lV +1 j 


LN(v + 1) 


- y = c 


LN(v + 1) 2-y + 2 c 2-y + 2 c 

#7: SOLVE y : c, v : (v : ■ - 1) v v = Y(e - 1)) 

l 2 J 

2 - y + 2-c 2-y + 2-c 

#8: DSOLVEl_GEHCy(e - 1), -1, x, y, k) = (ATAN(J(e - 1)) - x = -k) 

2-y + 2-c 2-y + 2-c 

#9: DS0LVE1_GEH(- - 1), - 1 , x, y, k) = (ATANG/Ce - 1)) + x = -k) 


2-y + 2-c 

#10: SOLVE(ATAMC f (e - 1» + x = -k, y) 

accepting some restrictions for the domain: 


#11: y = - 


LN(C0S(x + k) ) 


Example 21 is given in Fuster's paper. 
Example 16 is the last in my row: 
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2 ^ 

V 


#13: AUTONOMOUS 


v 


#14: DSOLVE1 


( V ) 

-GEN s -1, v s y = (LN(y) - LU(v) = -c) 

t v J 

u 


#15: SOLUTIONS (LN(y) - LN(v) = - Cj v > 


c c 

#16: DSOLVE1 _GEN (y ■ e 5 -1 , x s y s k) = (xe - LN(y) = k) 


#17: c 

SOLVE(x ■ e - LN(y) = L s y) = 


x ■ e - k 


ky = c 


) 


In my opinion Jose-Luis’ tools are a great support even in times of DERIVE 6. Referring to 
the original DERIVE tools requires a lot of more knowledge (which is not so bad) and of more 
manipulating. 

I wanted to compare DERIVE with the TIs and tried to solve Examples 9 through 24 using 
TI-NspireCAS (and the TI-92+ and Voyage 200): 


© Example 9 


deSolvely — ■y=x B sin(x),xj' 
lx 


y= -cos (x) -x's in (x)+ + c2 


© Example 10 


de So 1 ve (y 1 =>■' 1 vx+ iy" ' ) ^ , xjj 

-2 c5 (c5 I I 
3 m x m {e * I x | 1 J -cJj 


,e ’ ^sign(x)-e c -xj<0 and x>0 +c6 or y= 


-2 c? 


'x-(e c7 ‘|x|-l) c7 { . 


© Example 1 1 


de Solve (y'y-jl+yp ,xj;) 

ln(Jx 2 +2-c9\x+c9 2 -l +x+ c9j - (x+ c9) ■ Jx 2 + 2 ■ c9-x+ c9 2 - 


y= 


x 2 +2 - c9 - x+c 9 2 >1 and x+c9>0 +c10 c* 


© Example 12 

deSolvel 


H 

|-ln(cos(x+c17)),-cT1-^<x<yc1T+c72 


As there is the same CAS implemented it is no surprise that NSpire and the handhelds are 
behaving pretty the same (Examples 9, 10, 17 and 18): 


fri-TuriY FF? F ‘ |T Y rs Y ‘iT^i 

| t f — |fl 1 gebra |c a 1 c |other |Pr gn 1 0 |c 1 ean Up | | 

deSolue^y 1 1 --^-y 1 = x ■ sin(x) , x , yj 


y = "cos(x) - x ■ sin^x) + 


229 x 2 


+ 230 


1 deSolue(y 1 = y' 1 x+(y' ') 2 ,x,y) 

^ , -,- 2 e 33 . x .(^ 33 . |x| _ 

15 ■ xj 1 0 „ ^ - | + 23 ^ 


desolueCv 1 =w' IA 2,x^)l 

MHIH FiliD flUTD FUhK EE/30 


FF? YTF^t fht y F? t Fe^ t i 

| t f — |fl 1 gebra |c a 1 c |oihenPr gn 1 0 |c 1 ean Up | | 


1 deSolue(y 1 1 - (y 1 ) 3 + y (y 1 ) 3 = 0 , x ? y) 
y • ( y 2 - 3-y - 6-237) 


= x + 238 


1 deSolce(y 1 1 = (y 1 ) 3 ■ lnCy) ? x ? y) 


-y(2-y ln(y)-3-y + 4-239) 
^ = x + 240 


...solved 1 1 =w' 

MhIH FiliD HUTD FUhK EH/30 
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This is the rest of the list done with TI-NspireCAS: 
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Q 

y 2 = | Solve for v = y' 

1 + v 


Z.J 


c-y 

y 


v = i = ± ± 

dr \j v' 

dy-y i i t 

= = ax | Integrate again 

c-y 2 

s]c-y 2 =x + k | Solve for y 


I take the positive root and separate the variables 




y - ±-\Jc-(x + k ) 2 


I asked MATHEMATICA to do the job: 


in [2] := DSolve [y [x] * y 1 1 [x] + 1 + y * 1 [x] A 2 == 0 , y [x] , x] 

0ut[2]= -V -X 1 -2 c 2 - cl 

{m + }} 


Fortunately my result from above was confirmed. 


The next contribution causes problems. George Douros wrote his package for DERIVE 3 and 
it does not run with DERIVE 6 (solve - solutions and incompatibilities with the then used 
characters). I found George Douros email-address in the Internet and asked if there is a re- 
cent version available. I received an answer the following morning: 

Hello Josef Boehm! 

I am glad to hear from you... and sorry that I cannot be of help. 

I had written a much more powerful ODE package (at the same time with the Special Functions pack- 
age) for Derive 6. Some novelties in the simplification machine of Derive (I then had a few arguments 
with Albert Rich) prevented my new package from becoming functional. I have since (the end of 2003) 
stopped working on the package. 

Feel free to do whatever you think useful with the article or the package itself. My new demanding 
interests (ancient scripts...) do not allow me to make even a vague promise that I will be of some help 
to you. 

George Douros 


You can download the files and maybe that somebody can adapt the huge package - at 
least partially - for the recent DERIVE versions? 

I recommend the website demonstrating George's field of interest (ancient scripts). 


http://users.teilar.gr/~q1 951 d/ 
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Differential Equations with Derive ® 

George Douros 

Technological Education Institute of Larissa 
Larissa 41110, Greece 

Abstract 

This is a presentation of a new Derive package for solving Differential Equa- 
tions: ODE . MTH. The aim is to show how to use it and to explain some of the 
mathematics behind it, so that users can expand it, or use it as an educational aid. 

The package currently covers 1 st and 2 nd order Ordinary Differential Equations, 
both linear with variable coefficients and nonlinear. 

1. Introduction 

The first version of ODE . MTH appeared about a year ago and was written for Derive 2.5. Soft 
Warehouse found it extremely useful and encouraged the author to publicize and make it available to 
Derive users. In the mean time Derive 3 appeared, with new programming functions, and ODE .MTH 
grew into an entirely new package. Nonlinear equations are now covered almost exhaustively, singular 
solutions are found together with primitives, and a wide class of nonhomogeneous 2 nd order linear 
equations with variable coefficients that lead to solutions in terms of Bessel, Kummer and Hyper- 
geometric functions is treated in detail. 

The package, as well as this article, is based on Derive 3. The 4 User Manual , DERIVE Version 
3 ’ is the main reference. Other useful sources are listed in the Additional Resources appendix on pages 
331-340 of the Manual. 

2. The General Purpose Function: ODE (w, x, y , yl , y2 ) 

There are over 280 new functions defined in the package. Most of these are auxiliary, local 
functions. Some can be used independently. All are, however, auxiliaries to a single function, which 
solves ODEs of the 1 st and 2 nd order. This function, in terms of the default variables, is: 

ODE (w, X, Yr yl r y2 ) 

where, w is any equation (or function) involving the independent variable x, the unknown dependent 
variable y=y(x) and its derivatives yl=y'(x) and y2=y"(x). For example, to solve the equation 
y"(x) + y'(x)-3y(x)-2x -cosx, enter 

ODE (y2+yl-3*y=2*x A 2-COS (x) ) 


or, 


ODE (y2+yl-3*y-2*x A 2+COS (x) ) 
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To override the default variables in ODE all arguments must be entered explicitly. For exam- 
pie, to solve the equation u'(t ) + u{t) = t - 3 , enter 

ODE (yl+u=t A 2-3, t, u, yl) 


or, 

ODE (u_+u=t A 2— 3 r t, u, u_) 

Many of the methods in the package, use solutions of algebraic and/or transcendental equa- 
tions. It is, therefore, strongly recommended that the command: Manage-Branch-Any (this is actually 
the default DERIVE . INI setting) be issued before attempting to solve most ODEs which contain 
symbolic parameters. Exponents of symbolic powers x n of the independent variable x, appearing in 
the coefficients of y, y ' and y", in linear 2 nd order ODEs must be declared either positive or negative, 
since ODE needs to test the behaviour of the equation at its singular points. 

2.1. Solutions in terms of Special Functions 

Many 2 nd order ODEs will lead to solutions involving Special Mathematical Functions. These 
functions (Bessel, Kummer, Hypergeometric) are defined as arbitrary in the package and will not sim- 
plify. They can be simplified by loading the additional package FUN . MTH, after ODE . MTH. 


Loading both ODE and FUN in one session, will use up most of the available system memory 
in a plain 640K Derive session. In such cases the user can save the solutions, returned by ODE, and 
simplify them by loading FUN . MTH in a new session. No such problems occur with DeriveXM. 

2.2. Local Variables and Functions 

I know of no way to introduce local variables or functions in a Derive package. I have there- 
fore used its ability to ‘understand’ almost all ASCII characters in order to ‘imitate’ protected local 
variables and functions. This is done in the file ODE . ASC, where all variables are translated to ASCII 
characters above 180. The default variables (x,y,yl,y2) are left unchanged only in the ‘global 
functions’: ELIMINATE, RICCATI, GENERAL, ODE and the TRANS FORM_ utilities, to reduce typ- 
ing. If one uses only the default variables, it suffices, for example, to enter 

ODE ( (yl-y/x) A 2-x*y2) 

instead of 

ODE ( (yl-y/x) A 2-x*y2,x, y, yl, y2) 


3. Equations of the 1 st Order 

3.1. Equations linear in y f (x) 

The package uses the ‘standard’ techniques, but also looks for symmetries that allow simplify- 
ing transformations. For example, the equation = 0 is invariant under the 
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How did ODE proceed? When the solution algorithm reached the internal function CHANGE_X, 

M: 

equ_l was transformed by letting x^2arctan^ 2 (computed from the coefficients of y and y”), 

#19: INVERSE (INT ( V (DIF (equ_l , y) /DIF (equ_l , y2 ) ) ,x) ,x) 
x/ 2 

#20: 2 • ATAN ( e ) 

#21 : T RAN S F 0 RM_X (equ_l, #20) 

#22: 4* (y + y2) 

which is a simple equation to solve. 

4.2. Reducing an Equation to its Canonical Form 

The canonical form, y" +^-j(2qp f -2pq f + 4pr-q 2 jy = 0, of the differential equation 

equl:= py" + qy f + ry = 0 is obtained by letting y^ye^ , where $ -~\\^dx . 

It can easily be shown that equation equ2:= p{f y}" + q{f y}' + rf y = 0 has exactly the same 
canonical form as equ 1 and that its solution is a multiple of the solution of equ 1. 

This fact is used by ODE to solve equations like equ2 , by reducing them to their canonical 
form and proceed from that point by pattern matching with canonical forms of equations with known 
solutions. 

Equation hyp_ode is a simple case of a hypergeometric differential equation and is directly 
solved by ODE. 

#22 : hyp_: =2*x A 2* (1-x) *y2-x* (2*x+3) *yl+ (2*x+3) *y 
#23: ODE (hyp_) 

#24: y = x • ^cl_ • Vx • F2 1 , , , x j + c2_ • F2 1 ^2 , 0 , ' x ]] 

k y 

Suppose we transform hyp_ by first letting x — » x and then y —> . The resulting equa- 

Q{x) 

tion looks awesome. This is what we previously characterized as an unnatural form. By reducing this 
equation to its canonical form, however, ODE can solve it in a natural way, as if only by making sub- 
stitutions in #2 4. 

#25 : k A 2*TRANSFORM_X (hyp_, x A k) 

2 k k 2 k 

#26: 2*x • (1 - x ) * y2 - x*(2*x + 5*k - 2)*yl + k *(2*x + 3)*y 


#27 : Q (x) A 3*TRANSFORM_Y (#26, y/Q (x) ) 
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2 

k 

2x ( 1 

-x ) Q ( 


r k r 

xQ (x) 

2x 2 

r k r 

■ 2 

5 [ 

x Q (x) 

2 

rd n 

2x Q ( 

1 

X 

t: 

i 


2x — Q (x) -Q (x) -4x — Q (x) + (2-5k) Q (x) yl- 


Q (x) -2x 


\— Q(x)l 

'-dx -I 

!(x)j +X( 


2d 2 2- 

+xQ (x) — Q (x) + k Q (x) 
dx 


2 2 - 


Q(x)+4x — Q(x) +x(5k-2)Q(x) — Q(x)+3k Q(x) y 


Before solving the above equation, k must be declared positive (or negative) because ODE 
needs to test the behaviour of the equation at its singular points, and Derive is ‘reluctant’ to perform 
certain simplifications without this assumption for k. We also (optionally) declare x to be positive so 
that the solution #31 be fully simplified to look exactly like #2 4. 

#29: [x:eReal ( 0 , °°) , k : eReal (0,°°)] 

#30: ODE (#28) 

k r k/2 r513k-i r lk-i-i 


r k/2 r 5 

1 

3 k 

i r 1 k i i 

cl x F21 

r 

aX 

+ c2 F21 2,0, , x 

L 1-2 

2 

2 

J L 2 -l-l 


4.3. Searching for Symmetries 

The other major method used to bring an equation to its natural form is the search for symme- 
tries in the dependent and independent variables. Consider, for example, the readily integrable equa- 
tion #32, where Q(f) is an arbitrary function. If #32 is subjected to the transformations #3 3-# 3 7, 
the resulting equation #38 is in what we called an unnatural form. ODE can be impressive in solving 
it, only because it ‘realizes’ that #38 is invariant under the transformations (x,y) — » (<xx,/?y) and 
automatically reverses #3 3-# 3 7, to recover the original natural form. 

#32: y2-Q(yl) 

#33 : T RAN S F 0 RM_X (y2-Q (yl) , LN (x) ) 

2 

#34: - Q(yl-x) + y2 • x + yl-x 
#35 : T RAN S F 0 RM_ Y (#34, LN (y) ) 

r 2 

r yl‘ x i 2 y2 yl 

#36: - Q I + x * - 

L y J y 2 

L y 
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5. Additional Utilities 

The package makes internal use of some utilities that the end user may find useful. These utili- 
ties are: 

• ELIMINATE (v, a) tries to eliminate a from a parametric set v = [v l (a,...),v 2 (a , To sim- 
plify, for example, the solution of the equation xy' +4x = 2 yy\ which, as given by ODE, is: 

= = jc 1 (a 2 +2) ,y 2 - 4x 2 =0 , enter 

ELIMINATE ( [ [x=a*cl_, y=cl_* (a A 2+4) /2] , y A 2-4*x A 2=0] ,a) 

• TRANSFORM_X (w, t (x) , x, y, yl , y2 ) will transform the equation w, when x^Y(x). To 
transform, for example, the equation y'' + cot(x)y' + 4csc 2 (x)y = 0, by letting x^arccos(x), enter 

T RAN S F 0 RM_X (y2 + COT (x) *yl+4*CSC (x) A 2*y, ACOS (x) ) 

• TRANSFORM_Y (w, Q (x, y) , x, y, yl , y2 ) will transform the differential equation w, when 

y(x)—>Q(x,Y(x )) . To transform, for example, the equation y"-{x + l)y = 0, by letting 

1*2 

y^e 2 7(x), enter 

T RAN S F 0 RM_ Y (y2- (x A 2 + l ) *y, y*e A (x A 2/2) ) 

• TRANS FORM_P (w, x, y , yl , 0, r r rl ) will transform the 1 st order differential equation w, when 
(x,y(x))^r(#)-(cos#,sin<9). To transform, for example, the differential equation 

(y' 2 +1 )-(x-y) 2 =(x + yy') 2 , enter 

TRANS FORM_P ( (yl A 2+l) * (x-y) A 2= (x+y*yl) A 2) 

• RICCATI (w,s(x),x,y,yl) will return a solution of the Riccati DE w, which is more general 
than a known particular solution, s(x). For example, if y = ^ is a known solution of the Riccati e- 

quation y ' + \ + — = y 2 , to find a more general one, enter 

X x 

RICCATI (yl + l/x A 2+y/x=y A 2, 1/x) 
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• GENERAL (w, [1st (x) , 2nd (x) ] f x, y, yl, y2 ) will return the general solution of 
w — p(x)y" + q(x)y' + r(x)y + /(x) = 0, if one or two solutions of the homogeneous equation are 
known. If two solutions are known, they must be entered as a vector [ 1st (x) , 2nd (x) ] . If one 
solution is known, it may be entered as [lst(x),0],or simply 1st (x) . For example, to find 
the general solution of xy" (x) + xy' (x) = y(x) + x , if y=x is a known homogeneous solution, enter 

GENERAL (x*y2+x*yl=y+x, [x, 0] ) 


6. Conclusion 

ODE .MTH is more a mathematics than a programming package. I am a mathematician, not a 
programmer. I have tested it against similar packages in other symbolic mathematics programs and 
found it competitive. I will not make any further claims on its speed and efficiency and leave its eva- 
luation to those who will use it. This not simple modesty; it is the fear and dread of the easy counter- 
example to any of my claims. I have been faced with many of those in writing this package. With all 
due respect to the power given by the programming functions in the symbolic mathematics programs, I 
have been taught, once more, that a mathematical problem can be solved only by mathematics and not 
by programming functions. The latter are assistants; and Derive is A Mathematical Assistant par excel- 
lence! 


The hints, given here, on how ODE .MTH works are obviously not enough to help a user ex- 
pand, improve, or fully exploit it in a Differential Equations course, but they have made this ... a long 
article. 


ODE. MTH, together with its accompanying files, is available via the Derive Internet Mailing 
List. See the Additional Resources appendix of the Derive Manual. 


Office: 

Prof. George Douros 
T.E.I. of Larissa 
Larissa 41110 
Greece 


Home: 

George Douros 
Kolokotroni 3 
Larissa 41223 
Greece 


Tel: 041-611-061 ... 72 

Fax: 041-610-803 


Tel: 041-234-866 
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Asymptotically stable Solutions of the Systems 
of Ordinary Differential Equations 

A. Kozubik, Bratislava, Slovakia 

A given system of ordinary differential equations of the first order can be a mathematical 
model for a number of mechanic, biologic, economic etc. systems. The behaviour of this system is 
described by the solution of the respective system. For most of the systems we require that their be- 
haviour will be “similar” to the behaviour of some given system. This requirement is described exactly 
by the stability of the solution. 

Any system 

y'=g(t,y) (i) 

can be reduced by the transformation y = x + v where v is any solution of (1) into the system 

x' = f(t,x), f(t,o) = o (2) 

where 

f(t,x) = g(t,x + v)-g(t,v). 

System (2) has the trivial solution, which is the transform of the solution v of system (1). So, we can 
deal with the stability of the trivial solution of system (2). 

The trivial solution of system (2) is said to be stable iff for any x and any 8 there exists 
8 = 5(x,8) such that for all initial values S, satisfying || £, || < 8 and for all t > x the solution u(t;x£) of 
the initial problem 

x' = /(7,x), x(x) = £ 

satisfies the inequality 

|w(/;t, 0 I <e. 

The trivial solution of system (2) is said to be uniformly stable if number 8 does not depend on the 
initial point x. Moreover, the trivial solution is said to be uniformly asymptotically stable, if there ex- 
ists a number A > 0 such that for all || £, || < A the condition 

lim|| w(7;x,£,) || = 0 uniformly for all x 

holds. 

In this paper we deal with the linear system with constant coefficients in the form 

x f = Ax (3) 


where A is a real constant (w,^)-matrix. In this case we can apply the following assertions: 
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2-F(x2) 

x3 := x2 - , u [xl, x2 , x3, p + 1] 

ma 

#8: MULLER(mO, n) := APPEND([[xl, x2 , x3, p]], ITERATE 5 (u , m , mO, n)) 

3 

#9: F(x) := 1000 - x 

#10: MULLER([-0.5, 0, 0.5, 0], 12) 

xl x2 x3 p 

-0.5 0 0.5 0 

0 0.5 4000 1 

0.5 4000 0.80894382 2 

4000 0.80894382 1.1774703 3 

0.80894382 1.1774703 1.5253355 4 

1.1774703 1.5253355 17.43907 5 

# 11 : 

1.5253355 17.43907 8.2551095 6 

17.43907 8.2551095 9.6881509 7 

8.2551095 9.6881509 9.9871208 8 

9.6881509 9.9871208 10.000023 9 

9.9871208 10.000023 10 10 

10.000023 10 10 11 

10 10 10 12 

#12 : zl := 10 

3 

1000 - x 

#13: F(x) := 

x - zl 

#14: MULLERCE-0.5, 0, 0.5, 0], 5) 

xl x2 

-0.5 0 

0 0.5 

#15: 0.5 -5-8 .660254 - i 

-5 - 8 .660254 - i -5 - 8.660254 - i 
-5 - 8 .660254 - L -5 - 8.660254 - L 
. -5 - 8 .660254 - L -5 - 8.660254 - L 
#16: z2 := -5 - 8.660254- L 

3 

1000 - x 

#17: F(x) := 

(x - zl) - (x - z2) 


x3 p 

0.5 0 

-5 - 8.660254 - L 1 
-5-8 .660254 - i 2 
-5 - 8.660254 - i 3 
-5-8 .660254 - i 4 
-5-8 .660254 - i 5 
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#18: MULLER([-0.5, 0, 0.5, 0], 5) 




xl 


x2 



x3 

P 



-0.5 


0 



0.5 

0 



0 


0.5 

-5 

+ 

S. 660253999 - L 

1 

#19 : 


0.5 

-5 

+ 8. 660253999 - L 

-5 

+ 

S. 660254037 - L 

2 


-5 

+ 8. 660253999 ■ l 

-5 

+ S. 660254037 - l 

-5 

+ 

8.660254037 - l 

3 


-5 

+ S. 660254037 - i 

-5 

+ S. 660254037 - i 

-5 

+ 

8.660254037 - i 

4 


. -5 

+ S. 660254037 - L 

-5 

+ S. 660254037 - L 

-5 

+ 

8.660254037 - L 

5 


#20: [zl = 10, z2 = -5 - 8.660254- i. , z-3 = -5 + 8.660254- i.] 

3 

#21: F(x) := 1000 - x 


#22: MULLER ([-4, 2, 5, 0], 10) 


#23: 


xl 

x2 

x3 

P 

-4 

2 

5 

0 

2 

5 

15.85912687 

1 

5 

15.85912687 

9.266222135 

2 

15.85912687 

9.266222135 

9.934452212 

3 

9.266222135 

9.934452212 

9.999074829 

4 

9.934452212 

9.999074829 

10.00000014 

5 

9.999074829 

10.00000014 

10 

6 

10.00000014 

10 

10 

7 

10 

10 

10 

8 

10 

10 

10 

9 

10 

10 

10 

10 
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Some ideas for the revised Version of this DNL: 

(1) G P Speck writes about a “program” for demonstrating and performing Muller's Method. 
In 1995 it was a huge progress to have a list of functions - calling each other - working 
in the sense of a program. Now we can write real programs as a whole. I'll collect all 
procedures in one program without loosing the insight into the single steps. Application 
of ITERATES has an advantage in demonstrating the iterative nature of the process but 
the disadvantage that its syntax is sometimes difficult for students to understand. A sim- 
ple loop might be easier to follow. This realized in my program version. 
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nfi_meth(u, v, n, ca, cb, c, c 2 , ma, x4 , di’s, table) : = 

Pnog 

tables [ M xl M , M x2" , M x3 M , M p M ; v.j.1, v.j.2, v.j.3, 0] 

n_ := 1 

Loop 

If n_ > n 

RETURN table 

ca := (SUBST(u , x, vj,3) - SUBST(u, x, v4,2))/(v4,3 - v+2) 

cb := (SUBST(u, x, vj,2) - SUBST(u, x, vj,1))/(v 4,2 - vj,l) 

c := (ca - cb)/(v.j.3 - V 4 .I) 

c2 := ca + c- (vj, 3 - vj,2) 

dis := J(c2 A 2 - 4- SUBST(u, x, v.j.3). c) 

ma := IF(ABS(c2 + dis) > ABS(c2 - dis), c2 + dis, c2 - dis) 
x4 := v.j.3 - 2- SUBST(u, x, vj,3)/ma 
v := [v.j.2, v.j.3, x4] 


[12, 13, 14] 



table : = 
n_ :+ 1 

APPEND(table, 

[ [v J,1 , 

#2: 

(m_meth(1000 - 

3 

- x , [-0.5, 0 

, 0.5], 


" 9 .987121004 

10.00002340 

10 10 

#3: 

10.00002340 

10 

io n 


10 

10 

10 12 



3 


#4: r 

1000 - 
r n_meth 

- X 

, [-0.B, 

0, 0.B] 


-0.5 0 0.5 0 

0 0.5 -5 - S. 660254037 - i 1 

#5: 0.5 -5 - S. 660254037 - L -5 - S .660254037 - L 2 

-5 - S .660254037 - i -5 - S .660254037 - i -5 - S .660254037 - i 3 

-5 - S .660254037 ■ i -5 - S .660254037 - L -5 - S .660254037 - L 4 

. -5 - S .660254037 - i -5 - S .660254037 - i -5 - S .660254037 - i 5 

The next equation has only complex zeros. See the first steps 

4 2 

#10: m_meth(x + 7-x +10, [-10, 0, 10], 20) 


#11: FIRST (REVERSE (m_meth (x + 7-x +10, [-10, 0, 10], 20))) 

r -ios -ma -12s ] 

#12: |_7 .633765896-10 + 1.414213562- l , 4 . 034339544-10 + 1.414213562- i. , 2 .132092571-10 + 1.414213562- i , 20 J 





■"4 2 ^ 






x + 7-x +10 



#13: FIRST 

REVERSE 

m_nieth 

, [-10, 0, 10], 20 





V. 

x - 1.414213562- L J 

J 

J 


f -14 -14 -15 1 

#14: L - 5.066365831-10 - 1.414213561 - l , - 3.683475968-10 - 1.414213561 - l , - 1.0368848-10 - 1.414213561- l , 20 J 

Neglect the tiny real part, then the solutions are {-/V2, /V 2, -h/5, h/5}. 


In the next example the imaginary part can be taken as zero. 
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(m_meth(a - x , [-10, 0, 10], 20)) 


[ 20 , 21 , 22 ] 


-15 

-0. 70346742 + 1. 4981314-10 -i 

-0.70346742 

-19 

- 2.7887403-10 -i 

-0.70346742 

-28 

-6.5523512-10 -i 

IS 

-19 

-0. 70346742 - 2. 788 7403 -10 -L 

-0.70346742 

-28 

-6.5523512-10 -i 

-0.70346742 

-20 

- 2.0078468-10 -i 

15 

-28 

-0.70346742 - 6.5523512-10 -i 

-0.70346742 

-20 

- 2.0078468-10 -i 

-0.70346742 

-20 

- 3.9241406-10 -i 

20 


(m_meth(x - 10 + B - d. , [-10, 0, 10], 20)) 

[ 21 , 22 ] 

' 2,2094163 - 0,34420843 -l 2,2094163 - 0, 34420843 - L 2,2094163 - 0,34420843 -l 19 
. 2.2094163 - 0.34420843 - i 2.2094163 - 0. 34420843 - L 2.2094163 - 0. 34420843 - i 20 


(2) I wanted to illustrate the process by deriving the quadratic functions and plotting them 
together with the newly found approximation. 

2 

#24: q(x) := a-x + b-x + c 

3 

#25 : f (x) := 1000 - x 

Initial guesses are [-4,2,5] (see page 42): 

(5OLUTI0N5(q (-4) = f(-4) a q(2) = f(2) a q(5) = f (5) , [a, b, c])) 

#26: 1 

These are the coefficients of the parabola: 

#27: [-3, -18, 1040] 

You may use the FIT-function to obtain the parabola, too: 

( \ -4 f(-4) Ti 

r 2 I 2 

#28: FIT Lx, a-x + b-x + c 1 , 2 f(2) = - 3-x - 18-x + 1040 

, L 5 f(5) l 


#29: SOLUTI0NS(- 3-x - 18-x + 1040 =0, x) = [15.85912687, -21.85912687] 

Which zero of the approximating parabola is the appropriate one? 

#30: [ | 5 - 15.859126871, |5 - -21.85912687 | ] = [10.85912686, 26.85912686] 

We take the solution with the smaller distance to x3, which is 15.859... 

The next parabola should pass the points (2,f(2)), (5,f(5)) and (15.859...),f(15.859. 


#31: FIT 



2 

f (2) 


\ 2 1 




Lx, a-x + b-x + c I , 

5 

f(5) 



. 15.85912687 

f (15. 85912687) . 

J 


#32: - 22. 85912686 -x + 121. 0138880-x + 841.4087313 


#33: SOLUTIONSC- 22 ,85912686-x + 121. 013888-x + 841,4087313 , x) = [-3.972324047, 9,266222136] 


#34: [ | 15. 85912687 - -3. 972324047 | , | 15 ,85912687 - 9 . 266222136 | ] = [19.83145091, 6,592904733] 






D 


-N-L#20 


Additional Comments on Muller's Method 


p45 


This is the third parabola: 



r 2 i 

5 

f(5) 


#35: FIT 

Lx, a-x + b-x + c I , 

15 .85912687 

f (15. 8 591268 7) 



V 

. 9.266222136 

f (9. 266222136) . 

J 


2 

#36: - 30.12534900-x + 272 . BS09374-X + 265.2290376 

2 

#37: SOLUTIONSC- 30.125349-x + 272 . 5809374-x + 265.2290376, x) = [9.934452211, -0.SS62271639] 

The next figure shows the function together with the first three parabolas: 



(3) Finally I'd like to work with the Tl-Nspire using the spread sheet and the slider: 


/(x):=1000-x 3 Done | 

_y m 


1/99 



A 

B 

Cxi 

Dx2 

E X 3 

F ca 

G cb 

He 

■c2 

J dis 

Kma 

♦ 












1 

xl 

-0.5... 

-0.500... 

0 

0.500... 

0 

0 

0 

0 

0 

0 

2 

x2 

0 

0 

0.5000... 

4000.... 

-0.250... 

-0.25... 

0.000... 

-0.250... 

0.2500... 

-0.500... 

3 

x3 

0.50... 

0.5000... 

4000.0... 

0.808... 

-16002... 

-0.25... 

-4000... 

-3200... 

4472.3... 

-32006... 

4 



4000.0... 

0.8089... 

1.177... 

-16003... 

-1600... 

-4001 ... 

-1237.... 

41 86.6... 

-5424.... 

5 



0.8089... 

1.1774... 

1 .525... 

-2.993... 

-1600... 

-4001 ... 

-1477.... 

4262.1 ... 

-5739.... 

6 



1.1774... 

1.5253... 

17.43... 

-5.509... 

-2.99... 

-3.51 ... 

-6.730... 

118.50... 

-125.2... 

7 



1.5253... 

17.439... 

8.255... 

-333.0... 

-5.50... 

-20.1... 

-653.5... 

283.61 ... 

-937.1... 

8 



17.439... 

8.2551 ... 

9.688... 

-516.2... 

-333.... 

-27.2... 

-266.2... 

344.25... 

-610.5... 

9 



8.2551 ... 

9.6881 ... 

9.987... 

-241 .9... 

-516.... 

-35.3... 

-292.6... 

313.84... 

-606.5... 

10 



9.6881 ... 

9.9871 ... 

10.00... 

-290.3... 

-241.... 

-27.9... 

-298.7... 

299.43... 

-598.1... 

11 



9.9871 ... 

10.000... 

9.999... 

-299.6... 

-290.... 

-29.6... 

-299.9... 

299.99... 

-599.9... 

12 


9.9< 

399999996$ 

368*999... 

10.00... 

-300.0... 

-299.... 

-29.9... 

-299.9... 

299.99... 

-599.9... 

13 


_l 

|9.9999...| 

10.000... 

10.00... 

-300.0... 

-300.... 

-29.9... 

-300.0... 

300.00... 

- 600 . 0 ... y 


13 



9.9999... 

10.000... 

10.00... 

-300.0... 

-300.... 

-29.9... 

-300.0... 

300.00... 

-600.0... 

V 

14 



10.000... 

10.000... 

#UND... 

#UND... 

-300.... 

#UN... 

#UND... 

undef+... 

piece... 

15 



10.000... 

#UND... 

#UND... 

#UND... 

#UN... 

#UN... 

#UND... 

undef+... 

piece... 

16 



#UND... 

#UND... 

#UND„. 

#UND... 

#UN... 

#UN... 

#UND... 

undef+... 

piece... 

17 












C13 =d!2 
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The “undefined - messages” starting in row 14 are caused by the restricted accuracy. The 
values of x 2 and x 3 are so close that their difference is internally rounded to zero and then the 
division in column F gives an undefined result. All other error-messages are the conse- 
quence. 

However, we have result in the last complete row#13. 


I recommend the Document Settings as 
shown on the right: 


Document Settings () 


General General Settings 


Display Digits 

Fix 12 

V 

Angle 

Radian 

V 

Exponential Format 

Normal 

V 

Real or Complex Format 

Rectangular 

V 

Auto or Approximate 

Auto 


Vector Format 

Rectangular 

V 

Base 

Decimal 

V 

Unit System 

SI 

V 


Apply to 5ystem ][ Reset to Defaults ] [ OK ][ Cancel 


The entries for the cells are: 


Cl: =b1 

G2: 

=(f(d 1 )-f(c1 ))/(d 1 -cl ) 

CM 

_Q 

II 

Q 

H2: 

=(f2-g2)/(e1-c1) 

El: =b3 

12: 

=f2+h2*(e1-d1) 

C2: =d1 

J2: 

=V(i2 A 2-4*h2*f(e1 )) 

D2: =e1 

K2: 

when(abs(i2+j2) > abs(i2-j2),i2+j2,i2-j2) 

F2: = (f (e 1 )-f (d 1 ) )/(e 1 -d 1 ) 

E2: 

=e1-2*f(e1)/k2 


I try to perform the next step according to G P Speck's DERIVE routine in order to find the 
next (complex) root: 

I recommend deleting some rows in the table before changing the function and/or the initial 
guesses. The expressions in the cells might become too bulky and the system might hang 
up. (See below the complete entry of cell D6 - in grey.) 

It is better to start with less rows and then proceed by copying down row for row. Don't forget 
to save between the steps. 


0 

Done 


0 


1/2 


♦ 

A 

| 

n 

X 

°lx2 

ro 

X 

LU 

Bca 

Bcb 








1 

xl 

-0.5... 

-0.500000000000 

0 

0.500000000000 

0 


2 

x2 

0 

0 

0.500000000000 

-4.9999999999999-8.66025... 

-10.50... 

-9.5C 

3 

x3 

0,50... 

0.500000000000 

-4.9999999999999-8.6... 

-5.-8.6602540378443*7 

-5.500... 

-1 0.5 

4 



-4.9999999999999-8.... 

-5.-8.6602540378443*7 

-5.-8.6602540378443*7 

3.6602... 

-5.5C 

5 



-5.-8.6602540378443*7 

-5.-8.6602540378443*7 

-5.-8.6602540378443*i-(-1 .... 

#UND... 

3.66i 

6 



-5.-8.6602540378443*7 

-5.-8.6602540378443*... 

2*((3.4641016151377e-10-3... 

(59150... 

#UN 


■5.-3. 6602540378443*1- ( - 1 . 577350260 1 B96 e ■ 1 2 -4 . 226407308 1 038 e ■ 1 3*i)/(p i e c e wi s e (i*u n d ef+ u n d ef+ u n d e t u n d ef > u n d e t u n d ef- (i*u n d ef+ u 
ndefXundef<undef)) 


/fx): = 


lQQQ-x 

x-10 
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Finally I'd like to visualize the single steps of the algorithm by plotting the approximating pa- 
rabolas supported by a slider in the Graphs & Geometry Application. For this purpose I calcu- 
late the parabolas passing the three points using my polreg-program which is part of my per- 
sonal statistics library, called statistik (see DNL#...). 




DIx2 

EIx3 

®ca 

G cb 

H c 

1 c2 

J dis 

K rna 

L Iparabs 

M 

1 












30... 

2.000000... 

5.00000... 

0 

0 

0 

0 

0 

0 

-3.*x A 2-1 8.0Q0G0G0G0005*x... 


2 

30, 

5.000000... 

15.8591 ... 

-39.... 

-12.... 

-3.0... 

-48.... 

113.1... 

-16... 

-3 *x A 2-18. 000000000005*x+ 1040. 


3 

30, 

15.85912... 

9.26622... 

-35... 

-39.... 

-22.... 

-60... 

302.6... 

-90... 

-30.1 2534901 459*x A 2+ 272. 5... 


4 

12... 

9.266222... 

9.93445... 

-48... 

-35... 

-30.... 

-28... 

325.9... 

-61 ... 

-35.059801 21 9*x A 2+396. 560... 


5 

22... 

9.934452... 

9.99907... 

-27... 

-48... 

-35.... 

-30... 

304.5... 

-60... 

“29.1 9971 *x A 2+284.04306*x... 


6 

52... 

9.999074... 

10.0000... 

-29... 

-27... 

-29.... 

-29... 

299.9... 

-59... 

#ERR 



U 1 

=st atist ik\p 

'Olreg(2,{c 


4{f(< 



}) 






I introduce a slider and visualize the approximating parabolas step by step. This works pretty 
well until the accuracy of the system has reached its limit. 



/W:=e*-x 2 

Done 

1 1 



4/99 


k =3. 

—17— 

' y 

X 



B 

■xl 

Dlx2 

|A 

i 

* 





1 . 5 . 

/ 0.2 


i 

U, J... 

U.JUUUUUUUU... 

U = 

2 

2 

0 

0 

0.500000000... ■ 

3 

3 

0.50... 

0.5000000000... 

-0.717633299...- 

4 



-0.717633299... 

-0.702556274...- 

5 



-0.702556274... 

-0.703466489...- 


6.2 


6 



-0.703466489... 

-0.703467422...- 





7 



-0.703467422... 

-0.703467422...- 





8 







This is a 2 nd example (the equation mentioned by GP Speck on page 40). 


p48 


Thomas Weth: A Lexicon of Curves 


(7) 


D-N-L#20 


Ebene Algebraische und 
Transzendente Kurven (7) 

Thomas Weth, Wurzburg, Germany 


Die Kardiode - The Cardiod 

Legt man einen zylindrischen Ring auf eine ebene 
Tischflache, so werden einfallende Sonnenstrah- 
len am Inneren des Rings reflektiert. Die reflek- 
tierten Strahlen hiillen naherungsweise eine 
„Brennkurve“ ein, die wegen ihres Aussehens 
unter dem Namen Kardioide oder Herzkurve (vom 
griechischen icapdia) bekannt ist. Dieselbe Kurve 
lasst sich auch als helle Linie in einer mit Tee 
gefiillten Tasse beobachten. Genau genommen 
handelt es sich bei diesen Kurven um (halbe) 
Nephroiden. 

Ozanam erwahnt die Kardioide in seinem „ma- 
thematischen Worterbuch“ (Dictionnaire mathe- 
matique, Amsterdam, 1691); die heute noch giilti- 
ge Namensgebung geht allerdings auf Castillon 
zuriick (De curva cardioide, 1741 ). Die Kardioi- 
de ist mit der „fraktalen Geometrie“ zu neuer 
Aktualitat und Popularity gelangt - sie bildet (in 
grober Naherung) die Umrandung der "Ikone" der 
Chaostheoretiker, namlich der Mandelbrotmenge 
(volkstiimlich ,, Apfelmannchen“) . 




If you put a cylindrical ring on a plane table then 
the falling in sun rays will be reflected and they 
form approximately a "focal curve" which is called 
cardiod or heart curve . One can observe this curve 
in a cup filled with tee or coffee. Actually the car- 
diod consists of two (half) nephriods. In a rough 
approximation this curve forms the border of the 
chaos theorists' "icon" - the Mandelbrot set. 


Konstruktion der Kurve und Herleitung ihrer Gleichung 
Construction of the curve and derivation of its equation 



Eine genaue Konstruktion ergibt sich mit fol- 
gender Vorschrift: 

Konstruktion der Kardioide als katakaustische 
Linie - als Einhiillende einer Geradenschar. 
Gehen von einem Punkt auf dem Umfang eines 
Kreises Strahlen aus, die an der Kreislinie re- 
flektiert werden, so umhullen die reflektierten 
Strahlen die Kardioide. 

This is the instruction to find the cardiod as 

a catacaustic line: 

A point on the perimeter of a circle is the 
initial point of rays which are reflected at the 
perimeter. The family of the reflected rays 
form the cardiod as their envelope (see the 
sketch). 
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Now in 2009 we can produce the locus of the reflected segments on the graph screen of the 
TIs working with the Cabri-application. (Tl-Nspire offers the same occasion.) 






|L , . h 1 S r i 

|| 



[iff 

:J 1 IS 





The right graph shows the cardiod on the 
TI-NSpire's Graphs & Geometry Screen. 





WM 


l 

w 

MAIN DEG liUTD 

FUHl 



We proceed with Thomas Weth's text from 1995: 


Zunachst soil nun eine kartesische Koordinatendarstellung der Kurvenpunkte hergeleitet werden. 
Dabei ist zu beachten, dass die Kurvenpunkte genau die Beriihrpunkte mit den Kurventangenten sind. 
P sei der Punkt, von dem die Strahlen ausgehen. In obiger Figur wird der Strahl in Q reflektiert. Dann 

3a 

gilt fur die Steigung des Strahls QT: m = tan—. 


AuBerdem gilt nach dem Sinussatz im Dreieck OQT und unter Beriicksichtigung des Vorzeichens fur 

i 

also n = — 


den Achsenabschnitt n der Geraden QT: — = 

r 


. a 
sin — 
2 


sm 


n 3a 


. a 
r sin — 
2 _ 

3a * 


cos 


Fur die Kurventangente QT lautet also die Geradengleichung: y = x tan 


3a 

T 


. a 
r sin — 

o 

„ . Die partielle 

3a 
cos — 

2 


Ableitung nach dem Parameter fuhrt zur Gleichung der Hiillkurve. (Siehe Erklarung im Anhang.) 


The family of rays with their initial point on a fixed point on a circle form if they reflected at 
this circle the cardiod. Pay attention to the fact that the points of the curve are exactly the 
osculation points of the tangents. Let P the fixed point on the circle (the origin of the rays), 

Q is the intersection point of the ray and the circle. Then QT is the reflected ray with slope m. 
In triangle OQT we apply the sine rule to obtain the y-intercept n of the reflected ray. Using 
DERIVE we find the equation of the family of lines (I) with angle a as parameter. Partial dif- 
ferentiation wrt to the parameter (II) leads to the equation of the envelope. (See appendix 
and the accompanying file). 
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3a 

Formt man nun mit DERIVE um, so liefert die Multiplikation mit cos— : 


(I) 

(II) 


3a . 3a a 

Fix, v,a) = v cos x sin — + r sin — = 0. 

2 2 2 


df . 3a . 3a a 

-^— = -3 y sin 3x sin hr cos— = 0. 

da 2 2 2 


Lost man diese Gleichungen nach x undy auf, so erhalt man (nach mehreren Versuchen fiir die Ein- 
stellungen des Simplification Mode fiir Trigonometry und TrigPowers die Koordinaten der Kurven- 
punkte als: 

r r 

x = — (2 cos a -cos 2a) und y = — (2 sin a -sin 2a). 

Damit liefert DERIVE aus der Parameterdarstellung die gesuchte Kurve. 

2 2 2 
#1: x + y = r 

#2: P := [r-COS(t) , r-SIN(t)] 

#3: y = x-TAN 



#4: 


r r 

(2-C0S(t) - C0S(2-t)), (2- SIN(t) - SIN(2-t)) 

. 3 3 


# 5 : 


r r 

(2-C0SCO - C0S(2-O), (2-SINCO - SINC2-0) 

. 3 3 


#6: [[-r, 0], P] 

#7: [[0, 0], P] 



With DERIVE 6 we can improve the representation by introducing slider bars for the radius of 
the circle and for the angle which is formed by the ray PQ and the x-axis. 
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Zur Herleitung der Polar darstellung der Kurve bedarf es zunachst einiger Theorie: 

Gegeben sei eine Kurve (x(f), y(ff) mit dem Kurvenparameter t. Gesucht ist zu jedem 
Kurvenpunkt P der Abstand p zum Ursprung des Koordinatensystems und der Win- 
kel cp, den die Halbgerade OP mit der x-Achse einschlieBt (vgl. Skizze). 

Zu jedem ^benotigt man also zunachst den zugehorigen 
Parameterwert t. Erhalt man zwischen t und cp eine um- 
kehrbar eindeutige Beziehung (mit gewissen Differen- 
zierbarkeitseigenschaften), kann man t als Funktion von 
cp darstellen („t = t(cp)“) und kann die Kurve umparamet- 
risieren. Man erhalt dann (. x{t(cp )), y{t{cp))) oder kiirzer 
(x((p),y((p)). 

Prinzipiell kann man folgendermaBen vorgehen: 

Fiir den Winkel cp gilt: 

y(t ) t y(t ) 

<z> = arctan oder tan cp = . 

x(t) x{t) 



Diese Gleichung lost man - falls es gelingt - nach t auf und erhalt t in Abhangigkeit von cp, also eine 
Funktion cp — » t(cp). Einsetzen liefert dann (x(t((p)), y(t(<p))). Bildet man noch den Betrag p= \ (x(t(<p)), 
y{t{cp ))) |, so kann man mit DERIVE die Kurve in Polarform {p{cp),cp) zeichnen. 


Wendet man dieses Verfahren nun auf die oben hergeleitete Parameterdarstellung der Kardiode 

r r r \ 

— (2 cos a - cos 2a), — (2 sin a - sin 2a) an, erhalt man nach mehreren Umformungen eine Glei- 

y3 3 J 

chung 4. Grades mit der „Variablen“ cos (a), die DERIVE nicht losen kann. 


Verschiebt man allerdings die gesamte Kurve so, dass ihre Spitze in den Koordinatenursprung fallt, 

^ r r r \ 

— (2 cos a - cos 2a) — , — (2 sin a - sin 2a) und schlagt den oben angedeuteten Weg ein, erhalt 


3 


3 3 


J 


man: tan cp = 


sin<£> 
cos cp 


sm (x 

. Um den Radius fur die Polarform zu erhalten, kann man also bei dieser 

COStt 


speziellen Lage der Kurve direkt den Betrag der urspriinglichen Parameterdarstellung verwenden und 
erhalt (wieder nach einigen Umformungen mit DERIVE) die einfache Polardarstellung der Kardiode: 


(III) 


P = 


— (2 cos a - cos 2a) - — , —(2 sin a - sin 2a) 


2r 

= — (1-cos a). 
3 


Next we will find the polar form of this curve. We need for each point P(t) e curve its dis- 
tance p from the origin and the angle cp ( see the sketch above). If there exists a reversible 
unique relation between t and cp then t can be presented as a function of <p\ t = t(cp) which 

v(t) v(t) 

leads to {x{t((p)),y(t((p))) or(x(<p), y(<p)) with <p = arctan^-^ or tan<p = ^A 

x(t) x(t) 

If it is possible then we solve this equation for t and substitute t{cp) for t. With p = \(x(t(<p)), 
y{t((p)))\ we are able to plot this curve with DERIVE in polar form. 
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Calculating the length of the curve and the 
enclosed area gives interesting results. An 
interpretation of these results can be found 
investigating the cardiod as a special epi- 
cycloid. 

Merkwfirdige Ergebnisse ergeben sich beim Berechnen der Lange der Kurve und des von ihr einge- 
schlossenen Flacheninhalts. Mit dem CAS berechnet man: 




ist die Lange der Kurve ein ganzzahliges Vielfaches des zu 


ihrer Erzeugung verwendeten Kreises. 



p(oif da = 6r 2 7r , also wieder ein ganzzahliges Vielfaches der Flache des Kreises. Eine voll- 


standige Interpretation dieser Ergebnisse erhalten wir, wenn wir die Kardiode als eine spezielle Epi- 
zykloide, eine Kreisrollkurve betrachten. 


Anhang - Appendix 


Definition von Hfillkurven - Definition of envelopes 

Sei F(x,y,a) eine Funktion, so dass F(x,y,a) = 0 fur 
jedes a e ine Kurve beschreibt. F(x,y,a) bestimmt dann 
eine Kurvenschar mit dem Scharparameter a. (Im obi- 
gen Fall sind die F(x,y,a) die algebraischen Gleichun- 
gen von Geraden). Die Koordinaten der Schnittpunkte 
P zweier „benachbarter“ Kurven erffillen dann die Glei- 
chung 

F{x, y, a + Aa) - F(x , y,a) = 0 und speziell auch 

F Q, y,a + Aa)~ F(x, y,a) _ 

Aa 



dF (x y aA 

Fiir A a —> 0 geht die Gleichung fiber in: — — = 0. Erhalt man nun durch Elimination von a 

da 

dF (x y a ) 

aus den Gleichungen F(x,y,a) = 0 und — — = 0 wieder die Gleichung fix, y) = 0 einer Kurve, 

da 

so nennt man diese die Hiillkurve der Kurvenschar F(x,y,a). 


If F(x,y,a) is a family of curves with parameter a then we can find the envelope of these 
curves by eliminating the parameter from the two equations 

F(x,y,a ) = 0 and = o. If this procedure gives an equation f{x,y) of a curve then 

da 

this curve is the envelope. 
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I had the intention to place here a short Utility provided by Sergey Biryukov. As I have the strong 
impression that the contents of this issue has become “heavy enough”, I will do this in DNL#21. In- 
stead of this I reprint an international product: Belgian & Austrian. When I was in Belgium Jan Ver- 
meylen gave me a lot of work sheets and DERIVE files. Among them were some very interesting at- 
tempts to use the random number generator for creating exercise examples. These files were (be pa- 
tient, you will find them in one of the next year 's DNL)\ 

polyjrac.mth to practise factoring polynomials 

rekrij.mth to solve problems with arithmetic series 

quadr_eq.mth to create various quadratics 

At the DERIVE conference in Honolulu our working group had the idea to produce a “How it could 
look” work sheet for quadratics and factoring polynomials. We decided that I should add a demo file 
in combination with some exercises. I remembered Jan's quadr_eq.mth and this is the result: 

First load VIETAUTI . MTH as Utility file: 

#1: [zl := RAND0M(21) -10, z2 := RANDOM(31) - 15, z3 := RAN DOM (11) - 5] 

#2: list := [(a - zl)-(a - z2), (b - zl)-(b - z2), (d - zl)-(d - z2), (e - zl)-(e - 

z2) , (f - zl)-(f - z2) , (g - zl) • (g - z2), (h - zl).(h - z2), (i - zl).(i - z2), 
(k - zl) • (k - z2) , (o - zl) • (o - z2) , (p - zl).(p - z2), (t - zl).(t - z2), (x - 
zl) • (x - z2) , (y - zl) • (y - z2), (z - zl).(z - z2)] 

VIE(n, va) := [VECTOR (EXPAND( lim (x_ - zl).(x_ - z2)) =0, n_, 1, n)V 
L x_->va J 

VIETA(n) := EXPAND( [VECTOR (1 i st = 0, j, 1, n)l)' 

[ RAN DOM (15) + 1 J 

S01U(v) := VECTOR ( I F(DIM(S01UTI0NS(v_, (VARIABlES(v_)) )) = 2, S01UTI0NS(v_, 

1 

(VARIABlES(v_)) ), APPEND(S01UTI0NS(v_, (VARIABlES(v_)) ), S01UTI0NS(v_, 

1 1 

(VARIABlES(v_)) ))), v_ , v) 

1 

(Some details of the 1995 version for DERIVE 3 have been changed to make the file suitable for 
DERIVE 6 AO.) 

Then produce the following VIETADEM.DMO demo file. You can do it with any text editor or using 
the Edit > Annotation facility of the recent DERIVE versions. (More details can be found in DNL# 4.) 

The comments are the lines with the leading semi-colon. The must be entered as Annotation (without 
the semi-colon) together with the respective DERIVE commands: 

;This initializes the random number generator. Press ENTER. 

RANDOM ( 0 ) 

;Try to solve the equations. Continue with ENTER, 
task : =VIE (3,x) 

/Congratulations if you are right. Press ENTER for the next bundle. 

SOLU (task) 


#3: 

#4: 

#5: 






D 
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/Check your solutions pressing ENTER, 
task : -VIE (5, a) 

/Let's finish this sequence with a package of 10. 

SOLU (task) 

/ENTER will present all solutions, 
task : =VIETA ( 10 ) 

/If more than two answers are wrong, go on practising. You can then run the 
DEMO-file again. 

SOLU (task) 

/Interrupt the demo and check the solution using the SOLVE-command . 

VIETA ( 3 ) 

VIETA ( 5 ) 

/After checking the results, you may run the DEMO again. 

VIETA (10) 


If you are doing in with DERIVE 6 you might face problems saving the file as a “demo” file (with 
extensions .dmo only). I'd like to remind you that you have to save it as “vietadmo.dmo” - including 
the quotes!! 


2 

x - 17 -x = -72 


2 

g + g = 42 


VIETA(B) = 


2 

x + 7-x = S 


z = -7 v z = 7 


2 

. t - ll-t = -28 . 


Is it a feature or is it a bug? 

The function works - but if by chance - we have two solu- 
tions with z\ = -z 2 then DERIVE automatically solves the 
equation instead of only expanding the product of the two 
respective polynomials. 


DERIVE 5 & 6 allow programming. So we can collect the single procedures in a compact program 
which excludes the special cases z\ = - z 2 and z\ or z 2 = 0: 

#7: list_(zl, z2) :r [(a - zl)-(a - z2), (b - zl)-(b - z2), (d - zl)-(d - z2), (e - 

zl) • (e - z2), (f - zl) • (f - z2) , (g - zl).(g - z2), (h - zl).(h - z2), (i - 

zl).(i - z2) , (k - zl) • (k - z2) , (o - zl)-(o - z2), (p - zl).(p - z2), (t - 

zl).(t - z2) , (x - zl) • (x - z2) , (y - zl)*(y - z2), (z - zl)*(z - z2)] 


VIETA_(n, i_, tbl , p) := 

Prog 
i_ := 1 
tbl := [] 

Loop 

#8: If i_ > n 

RETURN tbl 
p := [zl, z2] 

If p ;1 + p;2 t 0 a pil*p|2 t 0 

tbl := APPEND(tbl , [ [EXPAND((1 ist_(pil, p *2)) i(RAND0M(15) + 1)) = 0]]) 
i_ :+ 1 
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Titbits from Algebra and Number Theory (6) 

by Johann Wiesenbauer, Vienna 


This time I would like to deal with a purely algebraic problem to make up for the fact that we have 
concerned mostly with number theoretic issues so far. And what could be more algebraic - at least 
from the historical point of view - than the solution of polynomial equations! As you all know 
DERIVE can easily cope with polynomial equations in one variable up to degree 4. Therefore to be a 
real challenge for DERIVE the degree of the polynomial in question has to be greater than 4. What 
about the polynomial equation 

z 17 -1 = 0 

which was solved by GauB when he was only 19 years old? (Actually GauB started his famous diary 
on March 30, 1796, with the following entry: “Principia quibus innititur section circuli, ac divisibili- 
tas eiusdem geometrica in septemdecim partes etc.” cf. [1]). Of course the application of the built-in 
SOLVE to this equation is definitely out of question. But what we could do is to follow the thoughts 
of genius GauB and enjoy their originality and elegance leaving all the drudgery to DERIVE. Are you 
ready for it? Then here we go. 

To begin with, it suffices to determine the solution 

z = e l<p = cos<z> + /sin<z> with <z>: — 

17 

of the equation above since all other solutions are merely powers of it. In the outline of the solution 
given by GauB in his letter to Gerling he actually determined cos cp which amounts to the same thing. 
To do so he first introduces the notations 


a:= cos (p + cos Acp 
Z?:=cos2<£> + cos8<£> 
c:=cos3<£> + cos5<£> 
d := cos6<£> + cos7<£> 

as well as e := a + b and f:=c + d. 

Then he states that “according to a wellkown theorem” the equality 


e + f = — 
2 


holds. Using 


cos k(p = ~[z k + z~ k ) = ~(z k + Z ll ~ k ) 


this can be derived (without DERIVE!) in the following way: 

8 1 16 i p 1 _i i 

e + f = Vcos kcp = —'S'z k = = — . 

h 2 tr 2 z-l 2 

GauB continues by claiming that 

ef=-l. 

Again we could calculate 

ef = (cos cp + cos 2 cp + cos 4 cp + cos 8 cp) (cos 3 cp + cos 5 cp + cos 6cp + cos lep) = 

= -(z + z l6 +z 2 +z l5 +z 4 +z n +z s +z 9 )-(z 3 +z u + z 5 +z u +z 6 + z n +z 7 + z 10 ' 
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and using z 17 = 1 as well as 




after a lot of tedious computations we could finally arrive at the desired result. (Try it out by yourself!) 
Is there a way to make DERIVE do these calculations for us? You won't be taken by surprise if I tell 
you there is actually one. The following utility function can be used quite generally to reduce a poly- 
nomial expression u in one variable modulo an equation v, which amounts to adding the rule v to the 
other rules for the handling of polynomials. 

RED(u , v) := ITERATE(RHS(v) ♦QUOTIENT(u_, LHS(v)) + REMAINDER(u_, LHS(v)) p u_ p u) 

2 4 8 9 13 15 16 

z + z + z +z + z + z +z +z 


By setting u := ef and taking (*) for v 
we get: 


5 6 7 10 11 12 14 

+ Z+ Z+ Z+ Z +z +z +z 


( 16 k_ A 

RED e *f p I z ~ -1 = -1 

l k_=l J 


Therefore e and / are both roots of the quadratic equation 


which has the solutions 


x 2 — (e + y*) x + e f-x 2 + — x — 1 — 0 


-l±Vl7 


Now we have come across a small problem, namely is 

e = x 1 ,f = x 2 (**) 

or vice versa? 

Again DERIVE can be of help in showing that the first alternative (**) is valid: 


APPROXf lim (COSO) + C0S(2-<p) + C0S(4*<p) + C0S(8*<p))) = 0.780776 
^ 2 - 71/17 

Do you see the achievement? The polynomial equation (*) of degree 16 has been broken up into two 
polynomial equations (**) with fewer terms! It takes no GauB to have the idea of trying the same trick 


once more: 


4 13 16 


+ z + z + 


2 8 9 15 

1 + z + z + z 


( 16 k_ A 1 

RED a*b p I z = -1 = 

L k_=l J 4 


Therefore a and b are both roots of the quadratic equation x 2 - {a + b)x + ab = x 2 - ex- — = 0. 
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Again we leave it to DERIVE to calculate its solutions: 

(2 1 T 

SOLUTIONS x - e-x = 0, x 

l 4 ) 


‘ 7(34 - 2 • Jl7 ) 7l 7 1 7(34 - 2«7l7) 7l 7 1 | 

+ , + 

8 8 8 8 8 8 . 

By applying APPROXIMATE one finds out that the first solution is a and the second one is b\ 

[1. 024740588, -0.2439641824] 

APPROX( lira (COS((p) + COSC4-*))) = 1.02474 
q>-*2 - jt/17 

In analogous manner the values of c and d are calculated: 

3 5 14 12 6 7 11 10 " 

z+z+z +z z+z+z +z 

c := , d := 

2 2 

r 16 k_ ^ 1 


RED c-d, 


f := - 


k_=l 

717 1 

4 4 


(2 1 

SOLVE x - f-x = 0, x 

L 4 


7C 34 - 2-V173 

x — — 

h 

’sj 

i 

M 

< 

1 

704 - 2-717) 

| 

717 

1 

8 

3 3 

3 

8 

8 

r 717 17 a 

X — f 1 

-/17 1 

f 717 17 ^ 

/ i 

7i? 

1 

L 32 32 J 

3 3 

■v i 

l 32 32 J 

3 

3 _ 


[x = 0.172075, x = -1.45235] 

APPROXC lim (CQS(3 -cfO + OOS(5-4))) = 0.172075 
4->2 -tt/17 

The rest is a piece of cake. Because of 


cos ^ cos V = ±(z + z 16 )i(z 4 + z 13 ) = i(z 5 + z 14 + z 20 + z 29 ) : 
= ^-(z 3 +z 14 +z 5 + z 12 ) = ^(cos3^ + cos5#?) = -^- 


we have 


(x - cos qj){x - cos 4 <p) = x 2 - (cos q> + cos 4<p)x + cos cp cos 4 <p - 

Q 

= x 2 - ax + — = 0. 

2 


We let DERIVE put the finishing touches: 
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a := -V 


17 J 17 


32 32 


/17 1 ( V17 17 

, c := J + 

8 8 32 32 


SOLUTIONS - a-x + =0, x 


s/C- 7(38-717 + 170) + 3-s/l7 + 17) /(34 - 2-^17) 


/(- /C38-/17 + 170) + 3-/17 + 17) /(34 - 2-^17) 


[0.9324722294, 0.09226835946] 

APPROXC Tim COS(<|>)) = 0.932472 
0-t2 -77/17 


That's it! There is only one thing that is slightly disturbing, namely that DERIVE didn't return the 
term for cos <p in its most beautiful form which is according to GauB: 

cosip = ~ + — ^Il7+— V34-2V17 +-Jl7 + 3Vl7-V34-2Vn -2^34 + 2^17. 

16 16 16 8 

But isn't this a bit of a tall order in view of this huge formula? 

Far more important is the fact that both forms of cos cp show that apart from the basic operations only 
square roots are used for its calculation. From this follows easily the stunning fact that the regular 
polygon with 17 edges can be constructed by using compass and straight-edge only! (cf [2] for the 
details!) As GauB was able to show this in his “Disqusitiones Arithmeticae” (1801) this is possible for 
a regular polygon with n edges if and only if n is of the form 

« = 2 ' . r>o,s>o, 

where fufi, . . .f s are pairwise distinct Fermat primes! Thus we’ve returned happily to number theory 
where we will continue next time! 


References 
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lagsgesellschaft, Leipzig 1981 
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(118-126) 

Three additional comments from the editor: 

Z> The file TITIBITS6.MTH is among the accompanying files and it contains the DERIVE 
functions and the complete calculation. 
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O Some time ago a DUG member sent a DERIVE file calculating the 17 roots of the equation. 
He refers to an article written by Friedrich Freytag, DdM 3, 1992 (pp 188 - 213). I recom- 
mend reading this article. There one can find a description how to construct the side of a 
regular 17-edge. I intend incuding this construction in DNL#21. 

O Johann Wiesenbauer has invented a remarkable function RED (u, v) . We compute a poly- 
nomial expression under consideration of an additional condition (which is in my opinion a 
kind of special substitution). Johannes told me that he had thought that this would be an 
advantage of other CAS packages, and he wanted to make this possible with DERIVE , too. 

Three easy to follow examples for applying RED (u, v) : 

[x := 2 - a + c , y := c + 3 -b] 


REDCx-y + 3 -x-y , a + b = 10) 

3 2 2 

- IS-b + 3 - b -(54 - c) + b - (4 - c + 123 -c + 130) + c-(c + 3)-(c + 20) 


[„= 


2 2 
3 - a + c,y:=c + 3-b 


RED 


3 - (x + y) , a - b = 


RED 


3 - (x + y) , a - b = 

2 J 


- 13-b + 9- b + 


3-Cfi-c + 3) 


RED(x + y, a - 2-b = 5) = 15-b + 60-b + 2-c + 75 

We would like to rewrite x + y = 3 a + 3 b 2 + 2c under consideration that a - 2b = 5. Let's reproduce 
Johann Wiesenbauer's RED-function step by step (= iteration by iteration): 

1 st step: Long division and its consequence 
(3a 2 + 3 b 2 + 2c ) : ( a - 2b) = 3a + 6b 

' 5 ' \3a 2 + 3 b 2 +2 c = 5(3 a + 6b) + I5b 2 +2c = 15 a + 15 b 2 + 30 b + 2c 

remainder : 15 b 2 + 2c 

2 nd step: We can perform the long division of the polynomials again: 

(15a + 15 b 2 + 30 b + 2c ) : (a - 2b) = 15 


Rem : 1 5 b 2 + 60 b + 2c 


\5a + \5b 2 + 30b + 2c = 3a 2 + 3b 2 +2c = l5-5 + \5b 2 + 60b + 2c 


No more division is possible, the next quotient would be zero and the result of the procedure remains 
the same 75 + 15 b 2 + 60 b + 2c. The iteration process has come to an end. Compare with the result of 
the RED-function from above! Things will change if we change the variable order: 

VariableOrder := [b, a] 


15 - a 


RED(x + y, a - 2-b = 5) 


15 - a S-c + 75 
+ 


(Who can write a respective function/program for the Tl's?) 


